Unitees : - ETM : mg/kg - Densitee apparentes : g/m3
Cette etude consiste a valoriser les donnees sur les Elements Traces Metalliques (ETM) issues de la premiere campagne (2000-2009) du Reseau de Mesure et de la Qualite des Sols (RMQS). Pour ce faire, plusieurs voies ont ete considerees, en commencant par une etude statistique sur les deux couches d’analyse du RMQS (H1 : 0-30 cm et H2 : 30-50cm) en utilisant les ETM selectionnes : Arsenic, Cadmium, Chrome, Cuivre, Cobalte, Molybdene, Nickel, Plomb et Zinc.
Les differentes pistes abordees sont les suivantes :
Dans le but d’obtenir une version plus recente de la table analyse volumetrique ponderee present dans l’entrepot (la derniere version datant de 2011), une extraction est effectuee sur la derniere version de la table mesures_echantillons_volumetriques du schema dm_donnees_ponctuelles. Cette table contient l’integralite des estimations massiques en elements grossiers (cailloux et graviers), en terre fine, le volume de l’echantillon…
Encore une fois, cette extraction est realisee uniquement sur les echantillons composites de la campagne 1. Deplus, les observations presentant des anomalies sont filtrees : Sites ne respectant pas la convention de prelevement par couche RMQS (prelevement par horizons de sol), observations avec une profondeur en sommet >50cm…
De plus, les echantillons preleves dans le cadre du projet BIOSOL n’ayant pas de donnees volumetriques pour la couche 2 (30-50cm), ne presente pas de stock ETM total. Neansmoins, un stock est calcule sur la couche 1 comme pour les autres points du RMQS.
Dans un soucis de simplification, et en attendant la mise a jour officielle des donnees ponderres de l’entrepot, les variables necessaires au calcul des stocks en ETM ont ete estimees comme suit :
\[Eg = (graviers+cailloux)/masse~seche\times100 \] Avec Eg le pourcentage en elements grossiers dans l’echantillon seche selon le rapport de la masse cumulee des fractions en elements grossiers presant dans l’echantillon.
# Commande SQL pour l'extraction des donnees volumetriques
connection=RODBC::odbcConnect("dela")
req_vol <- "SELECT id_campagne, id_site, type_profil_rmqs, no_horizon, no_prelevement,
prof_sup_moy, prof_inf_moy, prof_sommet, prof_base, id_analyse,
eg_graviers_masse as masse_eg_graviers,
eg_cailloux_masse as masse_eg_cailloux, da_volume_total as volume,
da_masse_humide as masse_humide, da_masse_seche as masse_seche,
da_no_methode, round((da_masse_seche/da_volume_total),2) as densite_apparente, round(( (eg_graviers_masse + eg_cailloux_masse) /da_masse_seche * 100 ), 2) as teneur_eg, round(((da_masse_humide-da_masse_seche)/da_masse_seche * 100), 2) as teneur_eau
FROM dm_donnees_ponctuelles.mesures_echantillons_volumetriques_rmqs
where id_campagne in ('1') and round((da_masse_seche/da_volume_total),2) is not null and type_profil_rmqs = 'C'
--and numero_annee_rmqs2 = '4'
--and region_partenaire ilike '%PACA%'
order by id_site, no_horizon, prof_sommet "
data_vol_tot <- sqlQuery(connection,req_vol)
odbcClose(connection) # fct de fermeture de la connection
# Nettoyage des donnees aberantes
Id_h3 <- unique(data_vol_tot[which(data_vol_tot$no_horizon%in%c(3,4,5,6,7,8)),"id_site"])
data_vol_tot_h3 <- data_vol_tot[which(data_vol_tot$id_site%in%Id_h3),]
data_vol_tot <- data_vol_tot[-which(data_vol_tot$id_site%in%Id_h3),]
data_vol_tot <- data_vol_tot[-which(data_vol_tot$prof_sommet>50),]
# Ponderation par id_site et par horizon (h1 - h2 - h9)
data_EG <- NULL
for(i in unique(data_vol_tot$id_site)){# traitement des donnees vol par site
dupli_ETM <- data_vol_tot[which(data_vol_tot$id_site==i),]
for(hh in unique(dupli_ETM$no_horizon)){# traitement des donnee par hz
dupli_ETM_hz <- dupli_ETM[which(dupli_ETM$no_horizon==hh),]
recap_meth <- as.data.frame(table(dupli_ETM_hz$da_no_methode))# Observation du nombre de methode dans l'horizon selectionne
if(nrow(recap_meth)>1){# si plusieurs methode appliquee
if(nrow(recap_meth[which(recap_meth$Var1==0.1),])>0){# On enleve les methode non identifie
recap_meth <- recap_meth[-which(recap_meth$Var1==0.1),]
}
if(hh==1 && c(2)%in%unique(dupli_ETM$no_horizon)){# Condition pour sauvegarder une continuite de la methode entre h1 et h2
recap_meth <- as.data.frame(table(dupli_ETM[which(dupli_ETM$no_horizon==2),]$da_no_methode))
}else if(hh==2){
recap_meth <- as.data.frame(table(data_EG[which(data_EG$no_horizon==1 & data_EG$id_site==i),]$no_methode))# respect de la continuite entre l'horizon h2 et l'horizon h1
}
if(length(unique(dupli_ETM$no_horizon))==1){
recap_meth <- recap_meth[1,]
}
dupli_ETM_hz <- dupli_ETM_hz[which(dupli_ETM_hz$da_no_methode==recap_meth[which(recap_meth$Freq==max(recap_meth$Freq)),1]),]
}
# Sauvegarde des epaisseur des echantillons volumetriques
Data_h_site <- as.data.frame(i)
colnames(Data_h_site) <- "id_site"
Data_h_site$no_horizon <-hh
Data_h_site$prof_sommet <- ifelse(is.na(sum(dupli_ETM_hz$prof_sommet)),min(na.omit(dupli_ETM_hz$prof_inf_moy)), min(na.omit(dupli_ETM_hz$prof_sommet)))
Data_h_site$prof_base <- ifelse(is.na(sum(dupli_ETM_hz$prof_base)),min(na.omit(dupli_ETM_hz$prof_sup_moy)), max(na.omit(dupli_ETM_hz$prof_base)))
Data_h_site$EG_pond <- (sum(na.omit(dupli_ETM_hz$masse_eg_graviers))+sum(na.omit(dupli_ETM_hz$masse_eg_cailloux)))/sum(na.omit(dupli_ETM_hz$masse_seche))*100
Data_h_site$DA_pond <- sum(na.omit(dupli_ETM_hz$masse_seche))/sum(na.omit(dupli_ETM_hz$volume))
Data_h_site$no_methode <- unique(dupli_ETM_hz$da_no_methode)
data_EG <- rbind(data_EG,Data_h_site)
}
}
data_ETM <- data_ETM[complete.cases(data_ETM$x_reel,data_ETM$y_reel),]
# Filtre les doublon de la table densite apparente
data_ETM$no_horizon <- as.numeric(as.character(data_ETM$no_horizon))
data_ETM <- left_join(data_ETM,data_EG,by=c("id_site"="id_site","no_horizon"="no_horizon"))
colnames(data_ETM)[87] <- "DA"
colnames(data_ETM)[86]<- "abondance_eg"
# Save data
write.csv2(data_vol_tot,paste0(Workdir,"/ETM_volumetrique_traitees.csv"))
write.csv2(data_ETM, paste0(Workdir,"/ETM_DA_pred_TOT.csv"),
row.names = FALSE)
# Importation du shapefile France
load("D:/WorkSpace/decoupage_administratif/france_entiere/bd_topo/contour_france.RData")
# Horizon 1
da_H1 <- data_ETM %>% filter(no_horizon == 1)%>%dplyr::select(abondance_eg)
mybreaks1 <- unique(round(unname(quantile(na.omit(da_H1$abondance_eg))),3))
ggplot() +
geom_polygon(data=FceMetCorse_coord, aes(long,lat, group = group), fill= "gray82", colour="grey") +
coord_equal() +
geom_point(data = data_ETM%>% filter(no_horizon == 1) ,
aes(x = x_reel, y = y_reel, size = abondance_eg,color=abondance_eg),shape=20, stroke=FALSE)+
scale_size_continuous(name="Element grossiers H1 (%)",breaks = mybreaks1,range=c(0.1,3.5))+
scale_color_viridis(name="Element grossiers H1 (%)",breaks = mybreaks1)+
theme_void() +
guides( colour = guide_legend()) +
theme(legend.position = "right")
# Horizon 2
da_H2 <- data_ETM %>% filter(no_horizon == 2)%>%dplyr::select(abondance_eg)
mybreaks2 <- unique(round(unname(quantile(na.omit(da_H2$abondance_eg))),3))
ggplot() +
geom_polygon(data=FceMetCorse_coord, aes(long,lat, group = group), fill= "gray82", colour="grey") +
coord_equal() +
geom_point(data = data_ETM %>% filter(no_horizon == 2),
aes(x = x_reel, y = y_reel, size = abondance_eg,color=abondance_eg),shape=20, stroke=FALSE)+
scale_size_continuous(name="Element grossiers H2 (%)",breaks = mybreaks2,range=c(0.1,3.5))+
scale_color_viridis(name="Element grossiers H2 (%)",breaks = mybreaks2)+
theme_void() +
guides( colour = guide_legend()) +
theme(legend.position = "right")
\[Da = massse~seche/volume~total \] Avec Da la densite apparente corespondant a la masse seche de l’echantillon par rapport a son volume total.
# Horizon 9
da_H9 <- data_ETM %>% filter(no_horizon == 9)%>%dplyr::select(DA)
mybreaks9 <- unique(round(unname(quantile(na.omit(da_H9$DA))),3))
ggplot() +
geom_polygon(data=FceMetCorse_coord, aes(long,lat, group = group), fill= "gray82", colour="grey") +
coord_equal() +
geom_point(data = data_ETM%>% filter(no_horizon == 9) ,
aes(x = x_reel, y = y_reel, size = DA,color=DA),shape=20, stroke=FALSE)+
scale_size_continuous(name="Densitee apparente H9 (g/cm3)",breaks = mybreaks9,range=c(0.1,3.5))+
scale_color_viridis( name="Densitee apparente H9 (g/cm3)",breaks = mybreaks9)+
theme_void() +
guides( colour = guide_legend()) +
theme(legend.position = "right")
# Horizon 1
da_H1 <- data_ETM %>% filter(no_horizon == 1)%>%dplyr::select(DA)
mybreaks1 <- unique(round(unname(quantile(na.omit(da_H1$DA))),3))
ggplot() +
geom_polygon(data=FceMetCorse_coord, aes(long,lat, group = group), fill= "gray82", colour="grey") +
coord_equal() +
geom_point(data = data_ETM%>% filter(no_horizon == 1) ,
aes(x = x_reel, y = y_reel, size = DA,color=DA),shape=20, stroke=FALSE)+
scale_size_continuous(name="Densitee apparente H1 (g/cm3)",breaks = mybreaks1,range=c(0.1,3.5))+
scale_color_viridis( name="Densitee apparente H1 (g/cm3)",breaks = mybreaks1)+
theme_void() +
guides( colour = guide_legend()) +
theme(legend.position = "right")
# Horizon 2
da_H2 <- data_ETM %>% filter(no_horizon == 2)%>%dplyr::select(DA)
mybreaks2 <-unique(round(unname(quantile(na.omit(da_H2$DA))),3))
ggplot() +
geom_polygon(data=FceMetCorse_coord, aes(long,lat, group = group), fill= "gray82", colour="grey") +
coord_equal() +
geom_point(data = data_ETM %>% filter(no_horizon == 2),
aes(x = x_reel, y = y_reel, size = DA,color=DA),shape=20, stroke=FALSE)+
scale_size_continuous(name="Densitee apparente H2 (g/cm3)",breaks = mybreaks2,range=c(0.1,3.5))+
scale_color_viridis(name="Densitee apparente H2 (g/cm3)",breaks = mybreaks2)+
theme_void() +
guides( colour = guide_legend()) +
theme(legend.position = "right")
Dans un soucis de transformation des unitees, nous utilisons la donnees DA, originellement en g/cm3, en kg/cm3. De plus, nous calculons l’epaisseur de chaque horizon a partir de la profondeur superieur et inferieur de l’horizon que nous gardons en cm.
# Preparation de la donnee
data_ETM$DA <- data_ETM$DA/1000 # g/cm3 en kg/cm3
data_ETM$Epaisseur <- abs(data_ETM$profondeur_hz_inf-data_ETM$profondeur_hz_sup) #calcul des epaisseurs en cm
# Function Recap
Stock_calc <- function(df,ETM,colout){
df[,colout] <- ifelse(is.na(df$abondance_eg),
df[,ETM]*df$Epaisseur*df$DA,# H9
df[,ETM]*df$Epaisseur*(1-(df$abondance_eg/100))*df$DA) # H1 et H2 mg/cm2
df[,colout] <- df[,colout]*10000 # mg/m2
return(df)
}
ETM.base <-c("as_tot_hf","cd_tot_hf","co_tot_hf","cr_tot_hf_RMQS","cr_tot_hf_RMQS1Bis","cu_tot_hf","ni_tot_hf",
"pb_tot_hf","zn_tot_hf","mo_tot_hf","cd_ext_edta","cr_ext_edta","ni_ext_edta","pb_ext_edta","zn_ext_edta",
"cu_ext_edta")
St.ETM <- c("Stock_as_hf","Stock_cd_hf","Stock_co_hf","Stock_cr_hf_RMQS","Stock_cr_hf_RMQS1Bis","Stock_cu_hf",
"Stock_ni_hf",
"Stock_pb_hf","Stock_zn_hf","Stock_mo_hf","Stock_cd_edta","Stock_cr_edta","Stock_ni_edta","Stock_pb_edta","Stock_zn_edta",
"Stock_cu_edta")
# Calcul des stock par horizon
for( j in 1:length(ETM.base)){
data_ETM <- Stock_calc(data_ETM,ETM.base[j],St.ETM[j])
}
write.csv2(data_ETM, paste0(Workdir,"/Stock_ETM_data.csv"),
row.names = FALSE)
Nous pouvons observer la variabilite de cette epaisseur pour chaque horizon :
# Epaisseur
mybreaks1 <- unique(unname(quantile(na.omit(data_ETM[which(data_ETM$no_horizon==1),"Epaisseur"])))[1:5])
mybreaks2 <- unique(unname(quantile(na.omit(data_ETM[which(data_ETM$no_horizon==2),"Epaisseur"])))[1:5])
mybreaks9 <- unique(unname(quantile(na.omit(data_ETM[which(data_ETM$no_horizon==9),"Epaisseur"])))[1:5])
# H1
ggplot() +
geom_polygon(data=FceMetCorse_coord, aes(long,lat, group = group), fill= "gray82", colour="grey") +
coord_equal() +
geom_point(data = data_ETM %>% filter(no_horizon == 1) ,
aes(x = x_reel, y = y_reel, size = Epaisseur,color=Epaisseur),shape=20, stroke=FALSE)+
scale_size_continuous(name="Epaisseur H1 (cm)",breaks = mybreaks1,range=c(0.1,3.5))+
scale_color_viridis(trans="log", name="Epaisseur H1 (cm)",breaks = mybreaks1)+
theme_void() +
guides( colour = guide_legend()) +
theme(legend.position = "right")
# H2
ggplot() +
geom_polygon(data=FceMetCorse_coord, aes(long,lat, group = group), fill= "gray82", colour="grey") +
coord_equal() +
geom_point(data = data_ETM %>% filter(no_horizon == 2) ,
aes(x = x_reel, y = y_reel, size = Epaisseur,color=Epaisseur),shape=20, stroke=FALSE)+
scale_size_continuous(name="Epaisseur H2 (cm)",breaks = mybreaks2,range=c(0.1,3.5))+
scale_color_viridis(trans="log", name="Epaisseur H2 (cm)",breaks = mybreaks2)+
theme_void() +
guides( colour = guide_legend()) +
theme(legend.position = "right")
# H9
ggplot() +
geom_polygon(data=FceMetCorse_coord, aes(long,lat, group = group), fill= "gray82", colour="grey") +
coord_equal() +
geom_point(data = data_ETM %>% filter(no_horizon == 9) ,
aes(x = x_reel, y = y_reel, size = Epaisseur,color=Epaisseur),shape=20, stroke=FALSE)+
scale_size_continuous(name="Epaisseur H9 (cm)",breaks = mybreaks9,range=c(0.1,3.5))+
scale_color_viridis(trans="log", name="Epaisseur H9 (cm)",breaks = mybreaks9)+
theme_void() +
guides( colour = guide_legend()) +
theme(legend.position = "right")
Afin de s’absoudre de l’heterogeneite des epaisseurs dans les donnees au sein des horizons, pour la spatialisation, nous avons choisis de calculer le stock total pour chaque ETM en cumulant les stocks de chaque couche (H1+H2+H9) :
\[\sum_{i=1}^n Stock_i = C_i \times Ep_i \times (1-Eg/100) \times Da_i\]
Avec C la teneur en ETM, Ep l’epaisseur de l’horizon, Eg le pourcentage d’elements grossiers et Da la densitee apparente.
Les cartogrammes presentes ci-dessous sont etablie a partir des coordonnees theoriques des points RMQS et ne correspondent pas a leur possition GPS reelle.
# CADMIUM
Stock_cd_HF <- data_ETM %>%
group_by(id_site,x_reel,y_reel,x_theo,y_theo,occupation,type_sol2,domaine_geo) %>%
dplyr::summarise(Sum_stock_cd_hf = sum(Stock_cd_hf))%>%
as.data.frame()
mybreaks1 <- unique(round(unname(quantile(na.omit(Stock_cd_HF[,"Sum_stock_cd_hf"]))),3))
test <- Stock_cd_HF%>%dplyr::select(x_theo,y_theo,Sum_stock_cd_hf)%>%
sf::st_as_sf(coords = c("x_theo","y_theo"), crs = 2154)%>%na.omit()%>%
mapview(zcol = "Sum_stock_cd_hf",layer.name = c("Stock Cadmium Total mg/m2"),cex="Sum_stock_cd_hf",at=mybreaks1,burst=FALSE, alpha = 0.5 )
Stock_cd_EDTA <- data_ETM %>%
group_by(id_site,x_reel,y_reel,x_theo,y_theo,occupation,type_sol2,domaine_geo) %>%
dplyr::summarise(Sum_stock_cd_edta = sum(Stock_cd_edta))%>%
as.data.frame()
save(Stock_cd_HF,file=paste0(StockDirectory,"/Stock_cd_HF.RData"))
save(Stock_cd_EDTA,file=paste0(StockDirectory,"/Stock_cd_EDTA.RData"))
# COBALTE
Stock_co_HF <- data_ETM %>%
group_by(id_site,x_reel,y_reel,x_theo,y_theo,occupation,type_sol2,domaine_geo) %>%
dplyr::summarise(Sum_stock_co_hf = sum(Stock_co_hf))%>%
as.data.frame()
mybreaks2 <- unique(round(unname(quantile(na.omit(Stock_co_HF[,"Sum_stock_co_hf"]))),3))
test2 <- Stock_co_HF%>%dplyr::select(x_theo,y_theo,Sum_stock_co_hf)%>%
sf::st_as_sf(coords = c("x_theo","y_theo"), crs = 2154)%>%na.omit()%>%
mapview(zcol = "Sum_stock_co_hf",layer.name = c("Stock Cobalte Total mg/m2"),cex="Sum_stock_co_hf",at=mybreaks2,burst=FALSE, alpha = 0.5 )
save(Stock_co_HF,file=paste0(StockDirectory,"/Stock_co_HF.RData"))
# CHROME
Stock_cr_edta <- data_ETM %>%
group_by(id_site,x_reel,y_reel,x_theo,y_theo,occupation,type_sol2,domaine_geo) %>%
dplyr::summarise(Sum_stock_cr_edta = sum(Stock_cr_edta))%>%
as.data.frame()
Stock_cr_hf_RMQS <- data_ETM %>%
group_by(id_site,x_reel,y_reel,x_theo,y_theo,occupation,type_sol2,domaine_geo) %>%
dplyr::summarise(Sum_stock_cr_hf_RMQS1 = sum(Stock_cr_hf_RMQS))%>%
as.data.frame()
mybreaks3 <- unique(round(unname(quantile(na.omit(Stock_cr_hf_RMQS[,"Sum_stock_cr_hf_RMQS1"]))),3))
test3 <- Stock_cr_hf_RMQS%>%dplyr::select(x_theo,y_theo,Sum_stock_cr_hf_RMQS1)%>%
sf::st_as_sf(coords = c("x_theo","y_theo"), crs = 2154)%>%na.omit()%>%
mapview(zcol = "Sum_stock_cr_hf_RMQS1",layer.name = c("Stock Chrome Total mg/m2"),cex="Sum_stock_cr_hf_RMQS1",at=mybreaks3,burst=FALSE, alpha = 0.5 )
save(Stock_cr_edta,file=paste0(StockDirectory,"/Stock_cr_edta.RData"))
save(Stock_cr_hf_RMQS,file=paste0(StockDirectory,"/Stock_cr_hf_RMQS.RData"))
# CHROME 1 et 1Bis
Stock_cr_hf_RMQS1<- data_ETM[which(data_ETM$no_horizon==1),
c("id_site","x_reel","y_reel","occupation","type_sol2","domaine_geo","Stock_cr_hf_RMQS")]
Stock_cr_hf_RMQS1Bis<- data_ETM[which(data_ETM$no_horizon==1),
c("id_site","x_reel","y_reel","occupation","type_sol2","domaine_geo","Stock_cr_hf_RMQS1Bis")]
save(Stock_cr_hf_RMQS1,file=paste0(StockDirectory,"/Stock_cr_hf_RMQS1.RData"))
save(Stock_cr_hf_RMQS1Bis,file=paste0(StockDirectory,"/Stock_cr_hf_RMQS1Bis.RData"))
# CUIVRE
Stock_cu_HF <- data_ETM %>%
group_by(id_site,x_reel,y_reel,x_theo,y_theo,occupation,type_sol2,domaine_geo) %>%
dplyr::summarise(Sum_stock_cu_hf = sum(Stock_cu_hf))%>%
as.data.frame()
mybreaks4 <- unique(round(unname(quantile(na.omit(Stock_cu_HF[,"Sum_stock_cu_hf"]))),3))
test4 <- Stock_cu_HF%>%dplyr::select(x_theo,y_theo,Sum_stock_cu_hf)%>%
sf::st_as_sf(coords = c("x_theo","y_theo"), crs = 2154)%>%na.omit()%>%
mapview(zcol = "Sum_stock_cu_hf",layer.name = c("Stock Cuivre Total mg/m2"),cex="Sum_stock_cu_hf",at=mybreaks4,burst=FALSE, alpha = 0.5 )
Stock_cu_EDTA <- data_ETM %>%
group_by(id_site,x_reel,y_reel,x_theo,y_theo,occupation,type_sol2,domaine_geo) %>%
dplyr:: summarise(Sum_stock_cu_edta = sum(Stock_cu_edta))%>%
as.data.frame()
save(Stock_cu_HF,file=paste0(StockDirectory,"/Stock_cu_HF.RData"))
save(Stock_cu_EDTA,file=paste0(StockDirectory,"/Stock_cu_EDTA.RData"))
# MOLYBDENE
Stock_mo_HF <- data_ETM %>%
group_by(id_site,x_reel,y_reel,x_theo,y_theo,occupation,type_sol2,domaine_geo) %>%
dplyr::summarise(Sum_stock_mo_hf = sum(Stock_mo_hf))%>%
as.data.frame()
mybreaks5 <- unique(round(unname(quantile(na.omit(Stock_mo_HF[,"Sum_stock_mo_hf"]))),3))
test5 <- Stock_mo_HF%>%dplyr::select(x_theo,y_theo,Sum_stock_mo_hf)%>%
sf::st_as_sf(coords = c("x_theo","y_theo"), crs = 2154)%>%na.omit()%>%
mapview(zcol = "Sum_stock_mo_hf",layer.name = c("Stock Molybdene Total mg/m2"),cex="Sum_stock_mo_hf",at=mybreaks5,burst=FALSE, alpha = 0.5 )
save(Stock_mo_HF,file=paste0(StockDirectory,"/Stock_mo_HF.RData"))
# PLOMB
Stock_pb_HF <- data_ETM %>%
group_by(id_site,x_reel,y_reel,x_theo,y_theo,occupation,type_sol2,domaine_geo) %>%
dplyr::summarise(Sum_stock_pb_hf = sum(Stock_pb_hf))%>%
as.data.frame()
mybreaks6 <- unique(round(unname(quantile(na.omit(Stock_pb_HF[,"Sum_stock_pb_hf"]))),3))
test6 <- Stock_pb_HF%>%dplyr::select(x_theo,y_theo,Sum_stock_pb_hf)%>%
sf::st_as_sf(coords = c("x_theo","y_theo"), crs = 2154)%>%na.omit()%>%
mapview(zcol = "Sum_stock_pb_hf",layer.name = c("Stock Plomb Total mg/m2"),cex="Sum_stock_pb_hf",at=mybreaks6,burst=FALSE, alpha = 0.5 )
Stock_pb_EDTA <- data_ETM %>%
group_by(id_site,x_reel,y_reel,x_theo,y_theo,occupation,type_sol2,domaine_geo) %>%
dplyr::summarise(Sum_stock_pb_edta = sum(Stock_pb_edta))%>%
as.data.frame()
save(Stock_pb_HF,file=paste0(StockDirectory,"/Stock_pb_HF.RData"))
save(Stock_pb_EDTA,file=paste0(StockDirectory,"/Stock_pb_EDTA.RData"))
# ZINC
Stock_zn_HF <- data_ETM %>%
group_by(id_site,x_reel,y_reel,x_theo,y_theo,occupation,type_sol2,domaine_geo) %>%
dplyr::summarise(Sum_stock_zn_hf = sum(Stock_zn_hf))%>%
as.data.frame()
mybreaks7 <- unique(round(unname(quantile(na.omit(Stock_zn_HF[,"Sum_stock_zn_hf"]))),3))
test7 <- Stock_zn_HF%>%dplyr::select(x_theo,y_theo,Sum_stock_zn_hf)%>%
sf::st_as_sf(coords = c("x_theo","y_theo"), crs = 2154)%>%na.omit()%>%
mapview(zcol = "Sum_stock_zn_hf",layer.name = c("Stock Zinc Total mg/m2"),cex="Sum_stock_zn_hf",at=mybreaks7,burst=FALSE, alpha = 0.5 )
Stock_zn_EDTA <- data_ETM %>%
group_by(id_site,x_reel,y_reel,x_theo,y_theo,occupation,type_sol2,domaine_geo) %>%
dplyr::summarise(Sum_stock_zn_edta = sum(Stock_zn_edta))%>%
as.data.frame()
save(Stock_zn_HF,file=paste0(StockDirectory,"/Stock_zn_HF.RData"))
save(Stock_zn_EDTA,file=paste0(StockDirectory,"/Stock_zn_EDTA.RData"))
test+test2+test3+test4+test5+test6+test7
# Importation des donnees sources
files.Rdata <- list.files(StockDirectory,pattern = paste0("Stock",".*",".RData$"))
for(i in 1:length(files.Rdata)){
setwd(StockDirectory)
Stock1 <- get(load(files.Rdata[i]))
if(i==1){
Stock_All <- Stock1
}else{
Stock_All <- left_join(Stock_All,Stock1[,c(1,ncol(Stock1))],by=c("id_site"))
}
}
Stock_All <- as.data.frame(Stock_All)
Stock_All <- Stock_All[-which(Stock_All$occupation=="foret"),]
# Creation d'un niveau de facteur de la covariable nom_occupation pour le cas de toutes occupations confondues pour le jeu de donnees dataSet
## Occupation
tmp <- Stock_All
tmp$occupation <- 'TtesOccupations'
Stock_All <- rbind.data.frame(Stock_All , tmp)
rm(tmp)
levels(Stock_All$occupation)[nlevels(Stock_All$occupation)] <- c("Toute_occupation_sauf_foret")
# Etudes statistiques relatives a la teneur en Elements chimiques ou en grandeurs physico-chimique
Grandeur <- c("Sum_stock_cd_edta","Sum_stock_cd_hf","Sum_stock_co_hf",
"Sum_stock_cr_edta","Sum_stock_cr_hf_RMQS1",
"Sum_stock_cu_edta","Sum_stock_cu_hf",
"Sum_stock_mo_hf","Sum_stock_pb_edta","Sum_stock_pb_hf","Sum_stock_zn_edta" ,"Sum_stock_zn_hf"
)
# Preparation des commentaires qui accompagneront les tables generees dans la partie suivante
GrandeurExplicite <- c("Stock total Cadmium EDTA","Stock total Cadmium HF","Stock total Cobalt HF","Stock total Chrome EDTA",
"Stock total Chrome RMQS1 HF","Stock total Cuivre EDTA","Stock total Cuivre HF",
"Stock total Molybdene HF","Stock total Plomb EDTA","Stock total Plomb HF","Stock total Zinc EDTA","Stock total Zinc HF")
# Definition de la variables de regroupement pour les tables finales
strateVar <- "occupation"
Stock_All$occupation <- as.factor(Stock_All$occupation)
levels(Stock_All$occupation) <- c("Culture permanente","Grande culture","Parc","Prairie","Sol nu","Toute occupation sauf foret","Vegetation naturelle","Zone humide")
strates <- levels(Stock_All[, strateVar])
Stock_All$no_horizon <- 1
couches <- 1
# Jointure du numero de cellule
Stock_All <- left_join(Stock_All,no_cellule[,c(9,10)],by=c("id_site"))
Stock_All <- Stock_All[-which(duplicated(Stock_All)),]
colnames(Stock_All)[ncol(Stock_All)] <- "no_cellule"
#setup parallel backend to use many processors
cores = detectCores()
cl <- makeCluster(3) #not to overload your computer
registerDoParallel(cl)
# #setup parallel backend to use many processors
# cores = detectCores()
# cl <- makeCluster(2) #not to overload your computer
# registerDoParallel(cl)
resufinal <-
foreach(
i = 1:length(Grandeur) ,#31:33 ,#
.combine = rbind.data.frame,
.export = "samplingDesignStats",
.packages = c("sp", "spcosa", "foreach"),
.verbose=FALSE
) %dopar% {
resucouche <-
foreach(
j = 1:length(couches),
.combine = rbind.data.frame
) %do% {
s=1
resParOccup <- foreach(strate = strates[c(1,2,4,6,7)],
.combine = rbind.data.frame
) %do% {
print(paste(
"Je traite ",
Grandeur[i],
"occupation",
strate,
" pour la couche" ,
couches[j],"----------------------------"
))
maskCouche <- Stock_All$no_horizon == couches[j]
maskStrate <- Stock_All[, strateVar] == strate
maskNA <- !is.na(Stock_All[, Grandeur[i]] )
if (sum(maskNA & maskCouche & maskStrate) > 30) {
resL <-
samplingDesignStats(grille,
Stock_All[maskNA & maskCouche & maskStrate , c("no_cellule", Grandeur[i]) ],
Grandeur[i])
# elt = Grandeur[i]
# dataSetSav = dataSet
# dataSet = dataSet[maskNA & maskCouche & maskStrate , c("no_cellule", Grandeur[i]) ]
resL <- as.data.frame(t(resL))
} else {
resL <- as.data.frame(t(rep(NA,15)))
colnames(resL) <- c("n_ind","Min","1stQuart","Median","Mean","3rdQu.","Max",
'1thDec','9thDec',
"Variance" , "se_mean_STSI",
"pcOutliers1.5","Vibrisse1.5","pcOutliers3","Vibrisse3")
}
cbind.data.frame(elt = GrandeurExplicite[i] ,couche = j, occup = strate , resL)
} # fin par occup
resParOccup
} # fin par couche
resucouche
} # fin par grandeur
stopCluster(cl)
# On sauve l'objet resufinal en cas de perte par ecrasementdes valeurs
SauvegardeResuFinal <- resufinal
# On ordonne les champs de la table resufinal
resufinal <- resufinal[ , c(1:5 , 11 , 6:9 , 12 , 10 , 13 , 14,15:18)]
# On renomme les champs de la table resufinal
colnames(resufinal) <- c("Grandeur_RMQS" ,
"Couche" ,
"Occupation" ,
"Effectif" ,
"Minimum" ,
"Premier_Decile" ,
"Premier_Quartile" ,
"Mediane" ,
"Moyenne" ,
"Troisieme_Quartile" ,
"Neuvieme_Decile" ,
"Maximum" ,
"Variance" ,
"Erreur_standard_STSI",
"pcOutliers1.5","Vibrisse1.5","pcOutliers3","Vibrisse3"
)
# Creation de la table xls
setwd(StockDirectory)
write.csv2(resufinal , "Stat_occup_Stock_tot_ETMRMQS18092020_.csv", row.names = FALSE)
kable(resufinal) %>%
add_header_above(c( "Statistique des stocks ETM (mg/m2) par occupation de sol" = 18))%>%
kable_styling(bootstrap_options = "striped", "bordered", full_width = F)%>%
scroll_box(width = "100%", height = "400px")
| Grandeur_RMQS | Couche | Occupation | Effectif | Minimum | Premier_Decile | Premier_Quartile | Mediane | Moyenne | Troisieme_Quartile | Neuvieme_Decile | Maximum | Variance | Erreur_standard_STSI | pcOutliers1.5 | Vibrisse1.5 | pcOutliers3 | Vibrisse3 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Stock total Cadmium EDTA | 1 | Culture permanente | 59 | 4.125722 | 15.551170 | 25.69607 | 41.76277 | 48.64609 | 59.70541 | 101.89190 | 139.7103 | 1.031030e+03 | 4.0300000 | 8.47 | 109.83890 | 0.00 | 139.7103 |
| Stock total Cadmium EDTA | 1 | Grande culture | 864 | 2.357347 | 21.900109 | 30.06454 | 45.96421 | 60.27719 | 72.53482 | 101.48264 | 1033.9086 | 5.125700e+03 | 2.2900000 | 4.40 | 134.46065 | 2.08 | 183.9732 |
| Stock total Cadmium EDTA | 1 | Prairie | 509 | 2.330283 | 18.039032 | 25.53477 | 39.61400 | 55.85873 | 63.39270 | 99.51503 | 1216.2817 | 5.911260e+03 | 3.3400000 | 6.48 | 119.73083 | 2.95 | 168.8127 |
| Stock total Cadmium EDTA | 1 | Toute occupation sauf foret | 1536 | 2.330283 | 18.647504 | 27.31468 | 42.27722 | 57.33511 | 67.85843 | 100.72304 | 1216.2817 | 5.082000e+03 | 1.7400000 | 5.01 | 127.03489 | 2.41 | 183.9732 |
| Stock total Cadmium EDTA | 1 | Vegetation naturelle | 90 | 2.964243 | 9.651386 | 13.25126 | 25.50497 | 45.77233 | 49.32581 | 96.19143 | 274.7434 | 3.020700e+03 | 5.4400000 | 10.00 | 93.47000 | 6.67 | 146.1918 |
| Stock total Cadmium HF | 1 | Culture permanente | 59 | 11.921656 | 31.603243 | 58.80555 | 94.63544 | 108.06390 | 143.62681 | 189.06954 | 313.7005 | 4.462680e+03 | 8.7600000 | 3.39 | 256.16789 | 0.00 | 313.7005 |
| Stock total Cadmium HF | 1 | Grande culture | 864 | 17.380621 | 43.212894 | 60.53393 | 98.76518 | 128.89107 | 156.18923 | 235.33344 | 1766.4990 | 1.622789e+04 | 3.9700000 | 4.86 | 296.26983 | 1.97 | 426.7281 |
| Stock total Cadmium HF | 1 | Prairie | 509 | 8.706180 | 35.869060 | 52.88167 | 78.11951 | 112.50206 | 123.39173 | 192.86506 | 3596.3162 | 3.481401e+04 | 8.1300000 | 6.48 | 226.18495 | 3.34 | 331.8238 |
| Stock total Cadmium HF | 1 | Toute occupation sauf foret | 1536 | 7.075067 | 37.793577 | 55.35645 | 88.16557 | 120.38569 | 146.70433 | 217.56010 | 3596.3162 | 2.159020e+04 | 3.5800000 | 5.34 | 281.76784 | 2.02 | 416.8914 |
| Stock total Cadmium HF | 1 | Vegetation naturelle | 90 | 7.940507 | 17.828598 | 31.91076 | 56.36471 | 97.42995 | 124.59676 | 233.58866 | 546.3982 | 1.074623e+04 | 10.1900000 | 6.67 | 253.30799 | 2.22 | 368.9678 |
| Stock total Cobalt HF | 1 | Culture permanente | 59 | 1000.677608 | 2249.318051 | 3149.42572 | 4850.19249 | 5296.56098 | 7037.68076 | 8932.18032 | 13885.7408 | 7.894492e+06 | 329.3700000 | 1.69 | 11494.44016 | 0.00 | 13885.7408 |
| Stock total Cobalt HF | 1 | Grande culture | 864 | 278.018333 | 1539.559871 | 2801.73206 | 4804.19883 | 5327.53649 | 6856.16395 | 9267.84317 | 44919.1643 | 1.463516e+07 | 122.9900000 | 3.70 | 12896.30245 | 0.81 | 18148.5041 |
| Stock total Cobalt HF | 1 | Prairie | 509 | 109.335513 | 1888.345856 | 3098.12992 | 5188.67839 | 6104.96541 | 7740.67600 | 10841.32099 | 37592.4344 | 2.052128e+07 | 199.4100000 | 3.73 | 14672.12473 | 1.38 | 21485.2750 |
| Stock total Cobalt HF | 1 | Toute occupation sauf foret | 1536 | 109.335513 | 1583.710225 | 2852.77634 | 4865.48823 | 5491.24113 | 7112.38759 | 9740.79752 | 44919.1643 | 1.609938e+07 | 88.6500000 | 3.84 | 13460.41938 | 1.04 | 19440.9276 |
| Stock total Cobalt HF | 1 | Vegetation naturelle | 90 | 397.025333 | 1226.363920 | 1890.48651 | 3459.46289 | 4033.28038 | 5405.79666 | 7716.21420 | 13883.6250 | 7.255172e+06 | 272.8800000 | 2.22 | 10656.02580 | 0.00 | 13883.6250 |
| Stock total Chrome EDTA | 1 | Culture permanente | 59 | 1.243404 | 14.484684 | 23.73079 | 39.39718 | 43.90060 | 61.00589 | 74.57792 | 143.9596 | 7.518600e+02 | 3.6500000 | 1.69 | 102.68994 | 0.00 | 143.9596 |
| Stock total Chrome EDTA | 1 | Grande culture | 864 | 2.823261 | 13.636152 | 21.08523 | 32.67109 | 39.10845 | 49.60732 | 69.96891 | 261.8277 | 7.611400e+02 | 0.8642523 | 4.75 | 92.22492 | 0.93 | 133.0533 |
| Stock total Chrome EDTA | 1 | Prairie | 509 | 2.479668 | 18.213482 | 28.42723 | 42.65084 | 53.63911 | 67.07076 | 92.57732 | 514.4019 | 2.190720e+03 | 2.1100000 | 3.93 | 124.92294 | 1.96 | 156.2955 |
| Stock total Chrome EDTA | 1 | Toute occupation sauf foret | 1536 | 1.243404 | 13.990087 | 22.99073 | 36.18522 | 44.59028 | 55.19771 | 82.37144 | 514.4019 | 1.359770e+03 | 0.8876911 | 4.95 | 103.42594 | 1.30 | 149.8895 |
| Stock total Chrome EDTA | 1 | Vegetation naturelle | 90 | 2.521620 | 7.793980 | 13.25745 | 27.68602 | 39.37031 | 56.42866 | 86.36474 | 249.1615 | 1.445150e+03 | 3.7900000 | 2.22 | 106.82381 | 1.11 | 153.7220 |
| Stock total Chrome RMQS1 HF | 1 | Culture permanente | 59 | 8409.209497 | 12189.047601 | 18234.14176 | 24522.75220 | 26362.68302 | 31776.83560 | 40710.00936 | 66558.1263 | 1.489958e+08 | 1572.1700000 | 3.39 | 49068.74040 | 0.00 | 66558.1263 |
| Stock total Chrome RMQS1 HF | 1 | Grande culture | 864 | 3208.285438 | 9593.594857 | 16394.32831 | 25797.00209 | 27893.67723 | 34046.70067 | 42482.51100 | 1181147.4495 | 1.761350e+09 | 1397.4200000 | 2.20 | 58959.40376 | 0.35 | 83813.5789 |
| Stock total Chrome RMQS1 HF | 1 | Prairie | 508 | 2681.345587 | 10879.095220 | 17791.22663 | 26729.21133 | 30663.53892 | 38205.70920 | 50076.68808 | 236463.1375 | 5.122649e+08 | 986.3700000 | 3.94 | 66649.39573 | 1.18 | 98355.8631 |
| Stock total Chrome RMQS1 HF | 1 | Toute occupation sauf foret | 1535 | 1628.037294 | 9742.986915 | 16391.95482 | 25578.43038 | 28554.13996 | 35183.12920 | 45201.54550 | 1181147.4495 | 1.267001e+09 | 879.6400000 | 3.06 | 63054.84840 | 0.91 | 87586.3689 |
| Stock total Chrome RMQS1 HF | 1 | Vegetation naturelle | 90 | 1910.435760 | 6349.600527 | 11270.40855 | 17853.19997 | 25717.14796 | 29824.61780 | 40294.40819 | 359486.7188 | 1.661218e+09 | 4096.3300000 | 4.44 | 46266.61514 | 2.22 | 64793.6731 |
| Stock total Cuivre EDTA | 1 | Culture permanente | 59 | 1025.808621 | 2032.176485 | 3559.55829 | 8582.68804 | 14927.86935 | 21173.11746 | 31282.33785 | 100304.8421 | 2.993281e+08 | 1972.3800000 | 5.08 | 42037.77667 | 1.69 | 60561.3102 |
| Stock total Cuivre EDTA | 1 | Grande culture | 864 | 112.123284 | 395.590294 | 652.46625 | 1167.87913 | 2175.62997 | 1877.97276 | 3723.75490 | 90416.7737 | 2.167034e+07 | 148.1200000 | 10.19 | 3656.21450 | 6.37 | 5504.0532 |
| Stock total Cuivre EDTA | 1 | Prairie | 509 | 42.868487 | 306.379880 | 559.23580 | 1032.36000 | 1849.71869 | 1868.97036 | 3249.93284 | 47039.4157 | 1.314052e+07 | 158.1600000 | 7.47 | 3826.74294 | 4.72 | 5796.3406 |
| Stock total Cuivre EDTA | 1 | Toute occupation sauf foret | 1536 | 42.868487 | 326.296260 | 602.23654 | 1119.84528 | 2545.22180 | 2002.20522 | 4565.96843 | 100304.8421 | 3.546391e+07 | 132.9800000 | 10.87 | 4099.69884 | 7.10 | 6122.5603 |
| Stock total Cuivre EDTA | 1 | Vegetation naturelle | 90 | 92.131346 | 165.481312 | 270.59940 | 623.62039 | 1794.50405 | 1134.71336 | 3299.57029 | 35344.6045 | 2.031673e+07 | 473.0800000 | 14.44 | 2164.87592 | 10.00 | 3154.8855 |
| Stock total Cuivre HF | 1 | Culture permanente | 59 | 5874.193263 | 8024.867168 | 12940.08650 | 24455.76590 | 36008.45229 | 53019.70730 | 69547.83796 | 152552.9327 | 9.736974e+08 | 3677.7200000 | 3.39 | 99697.27797 | 0.00 | 152552.9327 |
| Stock total Cuivre HF | 1 | Grande culture | 864 | 1013.736271 | 2893.510803 | 4598.67798 | 7052.80770 | 9201.96462 | 10635.20675 | 15131.35756 | 216289.4696 | 1.248605e+08 | 345.8000000 | 5.09 | 19628.91650 | 3.12 | 28273.8212 |
| Stock total Cuivre HF | 1 | Prairie | 509 | 435.992308 | 2850.974815 | 4831.02460 | 7682.95920 | 9147.11700 | 11589.45140 | 16386.92396 | 74772.5582 | 5.206143e+07 | 320.4700000 | 4.52 | 21611.99520 | 1.38 | 31088.3269 |
| Stock total Cuivre HF | 1 | Toute occupation sauf foret | 1536 | 318.716667 | 2791.942395 | 4667.11059 | 7399.34057 | 10208.92834 | 11417.66727 | 17489.86607 | 216289.4696 | 1.587101e+08 | 291.2300000 | 6.71 | 21272.43940 | 3.71 | 31599.6543 |
| Stock total Cuivre HF | 1 | Vegetation naturelle | 90 | 397.025333 | 1747.935708 | 3070.55668 | 5771.79774 | 8825.05160 | 10006.79342 | 18891.79978 | 67489.0900 | 1.070938e+08 | 1094.1300000 | 7.78 | 20555.86478 | 4.44 | 25610.9707 |
| Stock total Molybdene HF | 1 | Culture permanente | 59 | 90.414124 | 150.807595 | 227.64554 | 300.08452 | 414.08764 | 440.31419 | 683.26954 | 1907.2464 | 1.266157e+05 | 48.0100000 | 6.78 | 746.63022 | 6.78 | 746.6302 |
| Stock total Molybdene HF | 1 | Grande culture | 864 | 29.266599 | 112.980436 | 195.95981 | 281.41089 | 336.78792 | 386.87563 | 563.90055 | 3963.0574 | 7.893958e+04 | 9.1400000 | 6.25 | 670.70309 | 3.12 | 956.1433 |
| Stock total Molybdene HF | 1 | Prairie | 509 | 22.304445 | 144.615619 | 220.98736 | 315.56688 | 393.38614 | 450.41130 | 699.74039 | 3110.0726 | 1.052180e+05 | 14.1400000 | 7.66 | 793.32653 | 2.55 | 1135.5337 |
| Stock total Molybdene HF | 1 | Toute occupation sauf foret | 1536 | 22.304445 | 124.584557 | 202.88066 | 292.78622 | 363.56736 | 407.66050 | 611.26206 | 3963.0574 | 1.092007e+05 | 8.0400000 | 7.23 | 714.74314 | 3.91 | 1019.4917 |
| Stock total Molybdene HF | 1 | Vegetation naturelle | 90 | 24.890590 | 96.578952 | 149.77204 | 234.61261 | 434.12135 | 404.19044 | 809.88887 | 3923.9812 | 4.100332e+05 | 70.2500000 | 10.00 | 785.88252 | 6.67 | 1136.5337 |
| Stock total Plomb EDTA | 1 | Culture permanente | 59 | 390.977412 | 800.892506 | 1591.74009 | 2574.35947 | 3234.87824 | 4055.73484 | 6587.96388 | 10971.9977 | 5.424344e+06 | 302.8900000 | 6.78 | 7248.29583 | 0.00 | 10971.9977 |
| Stock total Plomb EDTA | 1 | Grande culture | 864 | 123.570628 | 874.962868 | 1364.18620 | 2012.29878 | 2693.21635 | 2982.41835 | 4317.64581 | 57143.1758 | 1.435340e+07 | 122.4200000 | 5.56 | 5408.50031 | 3.36 | 7616.9268 |
| Stock total Plomb EDTA | 1 | Prairie | 509 | 140.835291 | 858.411751 | 1465.62642 | 2347.84470 | 2907.94944 | 3457.25252 | 4784.06699 | 50770.7672 | 1.136706e+07 | 146.0000000 | 4.91 | 6395.15930 | 1.96 | 9364.7736 |
| Stock total Plomb EDTA | 1 | Toute occupation sauf foret | 1536 | 121.407456 | 835.619917 | 1367.22899 | 2097.93680 | 2791.41697 | 3214.32863 | 4622.41957 | 57143.1758 | 1.364693e+07 | 89.2600000 | 5.40 | 5985.76378 | 2.73 | 8679.8329 |
| Stock total Plomb EDTA | 1 | Vegetation naturelle | 90 | 121.407456 | 572.977264 | 1008.18015 | 1770.92280 | 2311.26906 | 2667.11082 | 4929.24528 | 15499.1747 | 4.957495e+06 | 221.3600000 | 8.89 | 5111.58440 | 2.22 | 7058.6638 |
| Stock total Plomb HF | 1 | Culture permanente | 59 | 2675.309303 | 6075.997345 | 11186.55141 | 16053.02360 | 16143.47000 | 20652.29255 | 23916.06335 | 37552.4264 | 5.337572e+07 | 929.8800000 | 3.39 | 28459.63480 | 0.00 | 37552.4264 |
| Stock total Plomb HF | 1 | Grande culture | 864 | 1201.486726 | 5870.608988 | 9014.80945 | 12432.91262 | 13851.93973 | 16394.39838 | 21847.69668 | 150073.9760 | 8.792232e+07 | 295.9000000 | 4.17 | 27438.38672 | 1.39 | 37294.0893 |
| Stock total Plomb HF | 1 | Prairie | 509 | 1895.877789 | 6465.066608 | 10332.64520 | 14279.55653 | 16587.90784 | 19027.13848 | 26649.64099 | 214497.8450 | 1.980525e+08 | 613.9300000 | 6.09 | 31583.05295 | 1.77 | 42803.6816 |
| Stock total Plomb HF | 1 | Toute occupation sauf foret | 1536 | 1201.486726 | 5891.588835 | 9129.78971 | 13038.78417 | 14757.71574 | 17264.52139 | 23815.19283 | 214497.8450 | 1.376464e+08 | 272.5400000 | 4.95 | 29462.13080 | 1.43 | 41120.8561 |
| Stock total Plomb HF | 1 | Vegetation naturelle | 90 | 1267.157303 | 3802.578913 | 6264.84448 | 8489.51600 | 10761.79734 | 12300.10586 | 19914.13849 | 47684.1523 | 6.164617e+07 | 785.5800000 | 8.89 | 19969.38489 | 3.33 | 28813.0188 |
| Stock total Zinc EDTA | 1 | Culture permanente | 59 | 113.775021 | 421.273870 | 623.16319 | 1136.95748 | 1646.83774 | 1818.69595 | 3369.51527 | 12996.3580 | 3.999035e+06 | 244.8700000 | 8.47 | 3579.66121 | 3.39 | 5326.2932 |
| Stock total Zinc EDTA | 1 | Grande culture | 864 | 98.978549 | 338.297157 | 489.87638 | 747.66753 | 1190.87379 | 1218.98879 | 2004.83390 | 94353.7178 | 1.214193e+07 | 115.5100000 | 7.75 | 2305.87576 | 2.78 | 3406.7599 |
| Stock total Zinc EDTA | 1 | Prairie | 509 | 87.673350 | 325.708172 | 489.56822 | 758.17160 | 1195.23247 | 1214.45471 | 2109.75653 | 48632.7083 | 6.554955e+06 | 112.9100000 | 8.45 | 2294.32372 | 3.93 | 2993.6444 |
| Stock total Zinc EDTA | 1 | Toute occupation sauf foret | 1536 | 87.673350 | 331.474461 | 493.18378 | 759.34148 | 1212.93693 | 1225.81809 | 2085.30210 | 94353.7178 | 9.400537e+06 | 77.1400000 | 8.27 | 2321.82433 | 3.65 | 3411.5323 |
| Stock total Zinc EDTA | 1 | Vegetation naturelle | 90 | 111.639060 | 268.829344 | 497.23469 | 717.11018 | 1114.06265 | 1128.81882 | 2082.57226 | 8801.9632 | 1.898698e+06 | 136.7100000 | 8.89 | 2085.73884 | 7.78 | 2393.9346 |
| Stock total Zinc HF | 1 | Culture permanente | 59 | 15319.087491 | 18227.737991 | 24082.38318 | 32436.48301 | 36942.81691 | 44607.82015 | 59489.54881 | 104208.1981 | 3.614133e+08 | 2472.2600000 | 5.08 | 71847.63333 | 0.00 | 104208.1981 |
| Stock total Zinc HF | 1 | Grande culture | 864 | 5134.264724 | 14656.482817 | 20853.42381 | 30363.58138 | 33799.43474 | 40932.87151 | 53466.31634 | 623538.2992 | 7.685910e+08 | 904.6900000 | 3.59 | 71036.97675 | 0.93 | 93355.3358 |
| Stock total Zinc HF | 1 | Prairie | 509 | 2440.726461 | 16418.620243 | 23917.60342 | 34613.74720 | 38989.14231 | 47365.56797 | 60977.60739 | 646923.1888 | 1.157353e+09 | 1507.5900000 | 3.14 | 80471.98917 | 0.59 | 111037.1620 |
| Stock total Zinc HF | 1 | Toute occupation sauf foret | 1536 | 1768.766667 | 15012.798333 | 21449.65629 | 31379.04293 | 35633.86026 | 43134.90885 | 57415.01668 | 646923.1888 | 9.127310e+08 | 761.1500000 | 3.78 | 75065.15200 | 0.85 | 104208.1981 |
| Stock total Zinc HF | 1 | Vegetation naturelle | 90 | 6684.732400 | 11767.364584 | 15232.40790 | 25272.00021 | 34212.29814 | 41706.42493 | 55740.42017 | 279695.2383 | 1.264904e+09 | 3701.8900000 | 4.44 | 76571.99200 | 2.22 | 87910.0186 |
# Importation des donnees sources
files.Rdata <- list.files(StockDirectory,pattern = paste0("Stock",".*",".RData$"))
for(i in 1:length(files.Rdata)){
setwd(StockDirectory)
Stock1 <- get(load(files.Rdata[i]))
if(i==1){
Stock_All <- Stock1
}else{
Stock_All <- left_join(Stock_All,Stock1[,c(1,ncol(Stock1))],by=c("id_site"))
}
}
Stock_All <- as.data.frame(Stock_All)
tmp <- Stock_All
tmp$domaine_geo <- 'TtesGeol'
Stock_All <- rbind.data.frame(Stock_All , tmp)
rm(tmp)
levels(Stock_All$domaine_geo)[nlevels(Stock_All$domaine_geo)] <- c("Toute_geol")
# Définition de la variables de regroupement pour les tables finales
strateVar <- "domaine_geo"
Stock_All$domaine_geo <- as.factor(Stock_All$domaine_geo)
levels(Stock_All$domaine_geo) <- c("Cristallin","Inconnu","Sedimentaire","Tourbiere","Tout domaine","Volcanique")
strates <- levels(Stock_All[, strateVar])
Stock_All$no_horizon <- 1
couches <- 1
# Jointure du numéro de cellule
Stock_All <- left_join(Stock_All,no_cellule[,c(9,10)],by=c("id_site"))
Stock_All <- Stock_All[-which(duplicated(Stock_All)),]
colnames(Stock_All)[ncol(Stock_All)] <- "no_cellule"
#setup parallel backend to use many processors
cores = detectCores()
cl <- makeCluster(3) #not to overload your computer
registerDoParallel(cl)
# #setup parallel backend to use many processors
# cores = detectCores()
# cl <- makeCluster(2) #not to overload your computer
# registerDoParallel(cl)
resufinal <-
foreach(
i = 1:length(Grandeur) ,#31:33 ,#
.combine = rbind.data.frame,
.export = "samplingDesignStats",
.packages = c("sp", "spcosa", "foreach"),
.verbose=FALSE
) %dopar% {
resucouche <-
foreach(
j = 1:length(couches),
.combine = rbind.data.frame
) %do% {
s=1
resParOccup <- foreach(strate = strates[c(1,2,3,5,6)],
.combine = rbind.data.frame
) %do% {
print(paste(
"Je traite ",
Grandeur[i],
"domaine_geo",
strate,
" pour la couche" ,
couches[j],"----------------------------"
))
maskCouche <- Stock_All$no_horizon == couches[j]
maskStrate <- Stock_All[, strateVar] == strate
maskNA <- !is.na(Stock_All[, Grandeur[i]] )
if (sum(maskNA & maskCouche & maskStrate) > 30) {
resL <-
samplingDesignStats(grille,
Stock_All[maskNA & maskCouche & maskStrate , c("no_cellule", Grandeur[i]) ],
Grandeur[i])
# elt = Grandeur[i]
# dataSetSav = dataSet
# dataSet = dataSet[maskNA & maskCouche & maskStrate , c("no_cellule", Grandeur[i]) ]
resL <- as.data.frame(t(resL))
} else {
resL <- as.data.frame(t(rep(NA,15)))
colnames(resL) <- c("n_ind","Min","1stQuart","Median","Mean","3rdQu.","Max",
'1thDec','9thDec',
"Variance" , "se_mean_STSI",
"pcOutliers1.5","Vibrisse1.5","pcOutliers3","Vibrisse3")
}
cbind.data.frame(elt = GrandeurExplicite[i] ,couche = couches[j], occup = strate , resL)
} # fin par occup
resParOccup
} # fin par couche
resucouche
} # fin par grandeur
stopCluster(cl)
# On sauve l'objet resufinal en cas de perte par écrasementdes valeurs
SauvegardeResuFinal <- resufinal
# On ordonne les champs de la table resufinal
resufinal <- resufinal[ , c(1:5 , 11 , 6:9 , 12 , 10 , 13 , 14,15:18)]
# On renomme les champs de la table resufinal
colnames(resufinal) <- c("Grandeur_RMQS" ,
"Couche" ,
"Domaine geol" ,
"Effectif" ,
"Minimum" ,
"Premier_Decile" ,
"Premier_Quartile" ,
"Mediane" ,
"Moyenne" ,
"Troisieme_Quartile" ,
"Neuvieme_Decile" ,
"Maximum" ,
"Variance" ,
"Erreur_standard_STSI",
"pcOutliers1.5","Vibrisse1.5","pcOutliers3","Vibrisse3"
)
# Creation de la table xls
resufinal <- na.omit(resufinal)
setwd(StockDirectory)
write.csv2(resufinal , "Stat_geol_Stock_tot_ETMRMQS18092020_.csv", row.names = FALSE)
kable(resufinal) %>%
add_header_above(c( "Statistique des stocks ETM (mg/m2) par domaine geologique" = 19))%>%
kable_styling(bootstrap_options = "striped", "bordered", full_width = F)%>%
scroll_box(width = "100%", height = "400px")
| Grandeur_RMQS | Couche | Domaine geol | Effectif | Minimum | Premier_Decile | Premier_Quartile | Mediane | Moyenne | Troisieme_Quartile | Neuvieme_Decile | Maximum | Variance | Erreur_standard_STSI | pcOutliers1.5 | Vibrisse1.5 | pcOutliers3 | Vibrisse3 | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 1 | Stock total Cadmium EDTA | 1 | Cristallin | 313 | 3.9989508 | 12.781617 | 17.95499 | 24.45503 | 27.59363 | 32.57603 | 44.52485 | 155.2103 | 2.645000e+02 | 8.122741e-01 | 4.79 | 54.08841 | 2.24 | 71.07806 |
| 2 | Stock total Cadmium EDTA | 1 | Inconnu | 32 | 5.4739805 | 8.101886 | 15.56010 | 27.35850 | 35.79023 | 40.32942 | 68.69351 | 146.3669 | 1.142650e+03 | 5.950000e+00 | 6.25 | 70.82124 | 6.25 | 70.82124 |
| 3 | Stock total Cadmium EDTA | 1 | Sedimentaire | 1309 | 2.3302833 | 21.607837 | 32.07823 | 49.99429 | 66.21804 | 76.44600 | 115.02100 | 1216.2817 | 6.084440e+03 | 2.080000e+00 | 5.81 | 142.99698 | 2.83 | 209.03843 |
| 4 | Stock total Cadmium EDTA | 1 | Tout domaine | 1683 | 2.3302833 | 16.959783 | 26.39146 | 41.97102 | 57.78189 | 68.68711 | 104.80991 | 1216.2817 | 5.054430e+03 | 1.630000e+00 | 5.76 | 129.73072 | 2.67 | 192.87670 |
| 6 | Stock total Cadmium HF | 1 | Cristallin | 313 | 4.4131278 | 26.469267 | 37.21533 | 50.86425 | 57.38985 | 68.27811 | 92.49688 | 288.1965 | 1.126540e+03 | 1.810000e+00 | 4.79 | 111.24720 | 1.92 | 135.38480 |
| 7 | Stock total Cadmium HF | 1 | Inconnu | 32 | 8.0222128 | 17.036503 | 23.27734 | 52.50760 | 67.20068 | 83.92248 | 132.23984 | 291.3307 | 4.080850e+03 | 4.630000e+00 | 6.25 | 133.73761 | 3.12 | 246.36137 |
| 8 | Stock total Cadmium HF | 1 | Sedimentaire | 1309 | 4.4518199 | 40.089391 | 62.75893 | 103.20555 | 136.20667 | 160.77204 | 247.25676 | 3596.3162 | 2.530592e+04 | 4.220000e+00 | 5.35 | 307.73334 | 2.52 | 450.38046 |
| 9 | Stock total Cadmium HF | 1 | Tout domaine | 1683 | 4.4131278 | 34.054113 | 52.91812 | 85.85400 | 119.14266 | 146.11333 | 221.98771 | 3596.3162 | 2.099939e+04 | 3.210000e+00 | 5.41 | 285.35808 | 2.14 | 420.71685 |
| 11 | Stock total Cobalt HF | 1 | Cristallin | 313 | 30.5634094 | 1339.932693 | 2204.74827 | 3902.07500 | 4665.93466 | 5962.02677 | 8285.74993 | 44919.1643 | 1.520392e+07 | 2.195600e+02 | 3.83 | 11531.77567 | 0.96 | 16303.32913 |
| 12 | Stock total Cobalt HF | 1 | Inconnu | 32 | 64.0198403 | 279.956853 | 704.66090 | 2502.86242 | 3578.10801 | 5025.25035 | 7787.98533 | 13883.6250 | 1.319755e+07 | 5.734500e+02 | 6.25 | 10525.14686 | 0.00 | 13883.62500 |
| 13 | Stock total Cobalt HF | 1 | Sedimentaire | 1309 | 6.3056939 | 1498.365034 | 2747.21468 | 4742.15893 | 5290.39851 | 7056.78611 | 9459.01792 | 33869.2661 | 1.349470e+07 | 9.006000e+01 | 3.06 | 13460.41938 | 0.84 | 18148.50411 |
| 14 | Stock total Cobalt HF | 1 | Tout domaine | 1683 | 6.3056939 | 1413.034291 | 2549.02326 | 4559.49722 | 5245.21844 | 6936.19803 | 9456.74788 | 44919.1643 | 1.580662e+07 | 9.130000e+01 | 3.57 | 13460.41938 | 0.89 | 20077.70494 |
| 16 | Stock total Chrome EDTA | 1 | Cristallin | 313 | 0.7426623 | 13.210388 | 24.16000 | 38.42818 | 46.12779 | 58.20287 | 80.04609 | 263.5351 | 1.228630e+03 | 1.660000e+00 | 4.79 | 102.93589 | 1.60 | 149.88946 |
| 17 | Stock total Chrome EDTA | 1 | Inconnu | 32 | 5.4490498 | 7.940649 | 21.29203 | 40.91753 | 52.55094 | 64.77952 | 92.59565 | 249.1615 | 2.315330e+03 | 8.840000e+00 | 3.12 | 133.06462 | 3.12 | 133.06462 |
| 18 | Stock total Chrome EDTA | 1 | Sedimentaire | 1309 | 0.6305694 | 14.801910 | 23.02846 | 36.78992 | 46.37584 | 56.79845 | 87.14055 | 514.4019 | 1.560100e+03 | 1.060000e+00 | 5.81 | 107.19878 | 1.68 | 157.03956 |
| 19 | Stock total Chrome EDTA | 1 | Tout domaine | 1683 | 0.6305694 | 14.016645 | 23.02388 | 36.84230 | 46.21482 | 56.97129 | 86.38377 | 514.4019 | 1.499250e+03 | 9.041181e-01 | 5.53 | 107.57127 | 1.54 | 158.39212 |
| 21 | Stock total Chrome RMQS1 HF | 1 | Cristallin | 313 | 147.1042602 | 7534.564672 | 12119.45164 | 20757.30066 | 25796.59620 | 29929.11888 | 37985.81623 | 1181147.4495 | 4.447744e+09 | 3.723470e+03 | 1.92 | 54814.20640 | 0.32 | 77079.87200 |
| 22 | Stock total Chrome RMQS1 HF | 1 | Inconnu | 32 | 1628.0372939 | 3268.123028 | 4407.87791 | 14454.10954 | 26291.69299 | 24739.58158 | 32389.36612 | 359486.7188 | 3.836724e+09 | 1.120414e+04 | 3.12 | 44534.77913 | 3.12 | 44534.77913 |
| 23 | Stock total Chrome RMQS1 HF | 1 | Sedimentaire | 1308 | 75.6683263 | 9230.490921 | 16290.95503 | 25233.26650 | 27174.49518 | 35177.96471 | 45755.80476 | 236463.1375 | 2.864512e+08 | 4.497100e+02 | 2.52 | 63054.84840 | 0.54 | 87586.36891 |
| 24 | Stock total Chrome RMQS1 HF | 1 | Tout domaine | 1682 | 75.6683263 | 8413.474957 | 14918.48796 | 24329.35881 | 27339.46789 | 34037.91931 | 44361.48896 | 1181147.4495 | 1.181885e+09 | 8.282700e+02 | 2.97 | 62450.02133 | 0.83 | 87586.36891 |
| 26 | Stock total Cuivre EDTA | 1 | Cristallin | 313 | 42.8684868 | 189.361654 | 372.50154 | 694.50360 | 1107.73789 | 1317.24539 | 2056.28928 | 12924.5355 | 2.290318e+06 | 8.542000e+01 | 5.43 | 2699.46230 | 3.51 | 3743.70289 |
| 27 | Stock total Cuivre EDTA | 1 | Inconnu | 32 | 81.7797985 | 149.071704 | 235.22684 | 725.05010 | 1614.53621 | 1555.46422 | 1955.52279 | 15804.7201 | 9.946262e+06 | 5.708700e+02 | 6.25 | 2850.80132 | 6.25 | 2850.80132 |
| 28 | Stock total Cuivre EDTA | 1 | Sedimentaire | 1309 | 27.4928252 | 324.812581 | 601.03078 | 1148.51619 | 2738.60643 | 2121.48932 | 4934.70510 | 100304.8421 | 4.081660e+07 | 1.535000e+02 | 11.61 | 4312.83804 | 7.56 | 6588.11217 |
| 29 | Stock total Cuivre EDTA | 1 | Tout domaine | 1683 | 27.4928252 | 285.357145 | 534.33271 | 1038.45865 | 2383.60299 | 1907.28455 | 4258.37274 | 100304.8421 | 3.281088e+07 | 1.244300e+02 | 10.40 | 3895.05584 | 6.83 | 6026.11420 |
| 31 | Stock total Cuivre HF | 1 | Cristallin | 313 | 164.2426206 | 2305.746970 | 3767.93395 | 6879.02693 | 8042.50820 | 11156.76800 | 15001.52342 | 31551.9708 | 3.089277e+07 | 2.946200e+02 | 2.88 | 20749.06275 | 0.00 | 31551.97082 |
| 32 | Stock total Cuivre HF | 1 | Inconnu | 32 | 463.2131292 | 730.955347 | 2175.44036 | 3918.07128 | 6689.30598 | 8455.77621 | 14599.80180 | 38004.7792 | 6.661682e+07 | 1.362200e+03 | 6.25 | 17199.28707 | 3.12 | 27944.28498 |
| 33 | Stock total Cuivre HF | 1 | Sedimentaire | 1309 | 71.7587961 | 2477.128437 | 4306.78549 | 7015.77241 | 10164.45354 | 11019.08650 | 17730.23910 | 216289.4696 | 1.820013e+08 | 3.323600e+02 | 7.33 | 20664.67421 | 4.51 | 30828.62434 |
| 34 | Stock total Cuivre HF | 1 | Tout domaine | 1683 | 71.7587961 | 2389.572942 | 4144.35490 | 6976.60886 | 9696.17037 | 11022.12266 | 17071.13634 | 216289.4696 | 1.499403e+08 | 2.687300e+02 | 6.36 | 21272.43940 | 3.51 | 31599.65426 |
| 36 | Stock total Molybdene HF | 1 | Cristallin | 313 | 5.7841966 | 102.394398 | 184.84891 | 267.63424 | 341.35913 | 372.59832 | 525.64708 | 3923.9812 | 1.421251e+05 | 2.100000e+01 | 6.71 | 617.55080 | 3.51 | 928.85800 |
| 37 | Stock total Molybdene HF | 1 | Inconnu | 32 | 37.6841250 | 57.734008 | 83.88618 | 175.90051 | 209.01522 | 282.54163 | 325.53741 | 832.4555 | 2.831602e+04 | 2.474000e+01 | 3.12 | 552.34482 | 0.00 | 832.45551 |
| 38 | Stock total Molybdene HF | 1 | Sedimentaire | 1309 | 5.0697779 | 114.173982 | 192.91963 | 288.72060 | 359.09760 | 408.70311 | 602.97425 | 4220.3441 | 1.083876e+05 | 8.590000e+00 | 6.49 | 725.90656 | 3.51 | 1053.76091 |
| 39 | Stock total Molybdene HF | 1 | Tout domaine | 1683 | 5.0697779 | 109.788138 | 189.99072 | 286.19639 | 355.41306 | 404.25690 | 600.09802 | 4220.3441 | 1.134374e+05 | 7.690000e+00 | 6.65 | 725.03371 | 3.39 | 1044.72660 |
| 41 | Stock total Plomb EDTA | 1 | Cristallin | 313 | 323.2880273 | 714.145930 | 1047.34336 | 1513.96026 | 1878.26702 | 2104.56119 | 3256.78869 | 15513.5624 | 2.470749e+06 | 8.393000e+01 | 6.71 | 3672.95280 | 3.19 | 5084.69568 |
| 42 | Stock total Plomb EDTA | 1 | Inconnu | 32 | 316.0448869 | 644.195843 | 1492.38066 | 2091.42253 | 3608.96297 | 2780.33436 | 4019.55456 | 46794.7066 | 6.427442e+07 | 1.386840e+03 | 9.38 | 4090.93137 | 6.25 | 5159.05683 |
| 43 | Stock total Plomb EDTA | 1 | Sedimentaire | 1309 | 89.3092402 | 965.074869 | 1562.20744 | 2364.57137 | 3041.95905 | 3452.49516 | 5102.86606 | 57143.1758 | 1.405538e+07 | 9.962000e+01 | 6.42 | 6284.17039 | 2.98 | 9058.13182 |
| 44 | Stock total Plomb EDTA | 1 | Tout domaine | 1683 | 89.3092402 | 847.170416 | 1375.87830 | 2104.66250 | 2798.43848 | 3227.91939 | 4684.50695 | 57143.1758 | 1.286525e+07 | 8.143000e+01 | 5.82 | 5985.76378 | 2.61 | 8766.02316 |
| 46 | Stock total Plomb HF | 1 | Cristallin | 313 | 575.5632703 | 6049.187975 | 8499.81953 | 12152.62560 | 13691.96634 | 16282.74390 | 21474.98049 | 74626.3200 | 7.626244e+07 | 4.791700e+02 | 4.47 | 27891.34210 | 2.24 | 36502.49590 |
| 47 | Stock total Plomb HF | 1 | Inconnu | 32 | 1710.6609375 | 3159.682985 | 5257.02236 | 9284.60926 | 15537.80103 | 15185.19926 | 22503.75951 | 161809.0198 | 7.693247e+08 | 4.700150e+03 | 6.25 | 28007.67763 | 3.12 | 32708.69833 |
| 48 | Stock total Plomb HF | 1 | Sedimentaire | 1309 | 563.7290308 | 5690.941866 | 8954.64819 | 12905.11650 | 14585.24297 | 17209.01700 | 23807.52078 | 214497.8450 | 1.309614e+08 | 2.980000e+02 | 4.66 | 29520.52787 | 1.22 | 41120.85615 |
| 49 | Stock total Plomb HF | 1 | Tout domaine | 1683 | 563.7290308 | 5588.903738 | 8558.29977 | 12541.24717 | 14312.61213 | 17039.87555 | 23287.87521 | 214497.8450 | 1.315790e+08 | 2.399600e+02 | 4.40 | 29593.91583 | 1.25 | 42320.26113 |
| 51 | Stock total Zinc EDTA | 1 | Cristallin | 313 | 87.6733500 | 291.681766 | 406.91592 | 595.55472 | 821.01171 | 952.03485 | 1516.30741 | 5374.0896 | 5.753318e+05 | 4.210000e+01 | 7.03 | 1740.74710 | 3.83 | 2501.92552 |
| 52 | Stock total Zinc EDTA | 1 | Inconnu | 32 | 155.2979186 | 339.657864 | 511.74989 | 679.56693 | 1333.35829 | 1116.03699 | 1828.06448 | 15132.2601 | 6.697846e+06 | 4.554800e+02 | 9.38 | 1834.15503 | 3.12 | 2838.36051 |
| 53 | Stock total Zinc EDTA | 1 | Sedimentaire | 1309 | 96.7487692 | 352.128658 | 521.33419 | 798.78291 | 1292.17580 | 1276.49309 | 2202.36392 | 94353.7178 | 1.085253e+07 | 8.820000e+01 | 8.63 | 2406.57603 | 3.97 | 3500.77901 |
| 54 | Stock total Zinc EDTA | 1 | Tout domaine | 1683 | 87.6733500 | 334.361826 | 493.98854 | 756.44071 | 1198.51815 | 1219.30413 | 2083.24789 | 94353.7178 | 8.710847e+06 | 6.997000e+01 | 8.26 | 2305.87576 | 3.74 | 3373.52440 |
| 56 | Stock total Zinc HF | 1 | Cristallin | 313 | 825.4976929 | 14819.217200 | 21520.34267 | 29929.51496 | 31755.67237 | 40931.11300 | 50577.33470 | 87839.0260 | 2.058327e+08 | 7.550200e+02 | 1.60 | 65677.17652 | 0.00 | 87839.02601 |
| 57 | Stock total Zinc HF | 1 | Inconnu | 32 | 2527.6722104 | 4699.826568 | 10781.52106 | 22612.65052 | 24183.04500 | 32248.66735 | 48003.65816 | 73760.2965 | 3.108414e+08 | 2.261410e+03 | 3.12 | 58327.53872 | 0.00 | 73760.29645 |
| 58 | Stock total Zinc HF | 1 | Sedimentaire | 1309 | 692.3651854 | 13034.881200 | 20216.90223 | 30786.76267 | 35338.18445 | 42886.16400 | 58532.17888 | 646923.1888 | 1.057430e+09 | 8.385800e+02 | 3.97 | 76571.99200 | 1.07 | 104208.19813 |
| 59 | Stock total Zinc HF | 1 | Tout domaine | 1683 | 692.3651854 | 13256.711183 | 20433.74352 | 30505.31020 | 34502.71504 | 42029.99226 | 56533.63622 | 646923.1888 | 8.762037e+08 | 6.943000e+02 | 3.86 | 74212.28733 | 0.83 | 104208.19813 |
# Importation des donnees sources
files.Rdata <- list.files(StockDirectory,pattern = paste0("Stock",".*",".RData$"))
for(i in 1:length(files.Rdata)){
setwd(StockDirectory)
Stock1 <- get(load(files.Rdata[i]))
if(i==1){
Stock_All <- Stock1
}else{
Stock_All <- left_join(Stock_All,Stock1[,c(1,ncol(Stock1))],by=c("id_site"))
}
}
Stock_All <- as.data.frame(Stock_All)
Stock_All <- dplyr::left_join(Stock_All,no_cellule[,c("id_site","code_dept")],by=c("id_site"))
Stock_All <- Stock_All[-which(duplicated(Stock_All)),]
Stock_All$Region <- ifelse(Stock_All$code_dept%in%c("29 ","22 ","56 ","35 "),"Bretagne",
ifelse(Stock_All$code_dept%in%c("50 ","14 ","76 ","27 ","61 "),"Normandie",
ifelse(Stock_All$code_dept%in%c("62 ","59 ","80 ","60 ","02 "),"Haut de France",
ifelse(Stock_All$code_dept%in%c("53 ","72 ","49 ","44 ","85 "),"Pays de la Loire",
ifelse(Stock_All$code_dept%in%c("95 ","78 ","91 ","77 ","75 ","93 ","92 ","94 "),"Ile de France",
ifelse(Stock_All$code_dept%in%c("28 ","45 ","41 ","37 ","36 ","18 "),"Centre val de Loire",
ifelse(Stock_All$code_dept%in%c("79 ","86 ","87 ","23 ","19 ","17 ","16 ","24 ","33 ","47 ","40 ","64 "),
"Nouvelle Aquitaine",
ifelse(Stock_All$code_dept%in%c("08 ","51 ","10 ","55 ","52 ","57 ","54 ","88 ","67 ","68 "),
"Grand Est",
ifelse(Stock_All$code_dept%in%c("89 ","58 ","21 ","71 ","70 ","39 ","25 ","90 "),
"Bourgogne Franche Comte",
ifelse(Stock_All$code_dept%in%c("03 ","63 ","15 ","43 ","42 ","69 ","07 ","01 ",
"26 ","38 ","74 ","73 "),"Auvergne Rhone Alpes",
ifelse(Stock_All$code_dept%in%c("46 ","12 ","48 ","30 ","34 ","81 ",
"82 ","31 ","32 ","65 ","09 ","11 ","66 "),
"Occitanie",
ifelse(Stock_All$code_dept%in%c("84 ","05 ","04 ","06 ","83 ","13 "),
"Provence Alpes Cote D'azur",
ifelse(Stock_All$code_dept%in%c("2A ","2B "),"Corse",NA)))))))))))))
tmp <- Stock_All
tmp$Region <- 'France entiere'
Stock_All <- rbind.data.frame(Stock_All , tmp)
rm(tmp)
levels(Stock_All$Region)[nlevels(Stock_All$Region)] <- c("France entiere")
# Définition de la variables de regroupement pour les tables finales
strateVar <- "Region"
Stock_All$Region <- as.factor(Stock_All$Region)
strates <- levels(Stock_All[, strateVar])
Stock_All$no_horizon <- 1
couches <- 1
# Jointure du numéro de cellule
Stock_All <- left_join(Stock_All,no_cellule[,c(9,10)],by=c("id_site"))
Stock_All <- Stock_All[-which(duplicated(Stock_All)),]
colnames(Stock_All)[ncol(Stock_All)] <- "no_cellule"
#setup parallel backend to use many processors
cores = detectCores()
cl <- makeCluster(3) #not to overload your computer
registerDoParallel(cl)
# #setup parallel backend to use many processors
# cores = detectCores()
# cl <- makeCluster(2) #not to overload your computer
# registerDoParallel(cl)
resufinal <-
foreach(
i = 1:length(Grandeur) ,#31:33 ,#
.combine = rbind.data.frame,
.export = "samplingDesignStats",
.packages = c("sp", "spcosa", "foreach"),
.verbose=FALSE
) %dopar% {
resucouche <-
foreach(
j = 1:length(couches),
.combine = rbind.data.frame
) %do% {
s=1
resParOccup <- foreach(strate = strates,
.combine = rbind.data.frame
) %do% {
print(paste(
"Je traite ",
Grandeur[i],
"domaine_geo",
strate,
" pour la couche" ,
couches[j],"----------------------------"
))
maskCouche <- Stock_All$no_horizon == couches[j]
maskStrate <- Stock_All[, strateVar] == strate
maskNA <- !is.na(Stock_All[, Grandeur[i]] )
if (sum(maskNA & maskCouche & maskStrate) > 30) {
resL <-
samplingDesignStats(grille,
Stock_All[maskNA & maskCouche & maskStrate , c("no_cellule", Grandeur[i]) ],
Grandeur[i])
# elt = Grandeur[i]
# dataSetSav = dataSet
# dataSet = dataSet[maskNA & maskCouche & maskStrate , c("no_cellule", Grandeur[i]) ]
resL <- as.data.frame(t(resL))
} else {
resL <- as.data.frame(t(rep(NA,15)))
colnames(resL) <- c("n_ind","Min","1stQuart","Median","Mean","3rdQu.","Max",
'1thDec','9thDec',
"Variance" , "se_mean_STSI",
"pcOutliers1.5","Vibrisse1.5","pcOutliers3","Vibrisse3")
}
cbind.data.frame(elt = GrandeurExplicite[i] ,couche = couches[j], occup = strate , resL)
} # fin par occup
resParOccup
} # fin par couche
resucouche
} # fin par grandeur
stopCluster(cl)
# On sauve l'objet resufinal en cas de perte par écrasementdes valeurs
SauvegardeResuFinal <- resufinal
# On ordonne les champs de la table resufinal
resufinal <- resufinal[ , c(1:5 , 11 , 6:9 , 12 , 10 , 13 , 14,15:18)]
# On renomme les champs de la table resufinal
colnames(resufinal) <- c("Grandeur_RMQS" ,
"Couche" ,
"Region" ,
"Effectif" ,
"Minimum" ,
"Premier_Decile" ,
"Premier_Quartile" ,
"Mediane" ,
"Moyenne" ,
"Troisieme_Quartile" ,
"Neuvieme_Decile" ,
"Maximum" ,
"Variance" ,
"Erreur_standard_STSI",
"pcOutliers1.5","Vibrisse1.5","pcOutliers3","Vibrisse3"
)
# Creation de la table xls
resufinal <- na.omit(resufinal)
setwd(StockDirectory)
write.csv2(resufinal , "Stat_region_Stock_tot_ETMRMQS18092020_.csv", row.names = FALSE)
kable(resufinal) %>%
add_header_above(c( "Statistique des stocks ETM (mg/m2) par region" = 19))%>%
kable_styling(bootstrap_options = "striped", "bordered", full_width = F)%>%
scroll_box(width = "100%", height = "400px")
| Grandeur_RMQS | Couche | Region | Effectif | Minimum | Premier_Decile | Premier_Quartile | Mediane | Moyenne | Troisieme_Quartile | Neuvieme_Decile | Maximum | Variance | Erreur_standard_STSI | pcOutliers1.5 | Vibrisse1.5 | pcOutliers3 | Vibrisse3 | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 1 | Stock total Cadmium EDTA | 1 | Auvergne Rhone Alpes | 222 | 5.550854e+00 | 15.418083 | 22.11181 | 32.85589 | 44.92847 | 53.31763 | 88.15434 | 227.8665 | 1.460820e+03 | 2.280000 | 7.66 | 97.43808 | 3.15 | 146.19180 |
| 2 | Stock total Cadmium EDTA | 1 | Bourgogne Franche Comte | 133 | 5.337040e+00 | 23.704711 | 37.33545 | 55.98297 | 84.34368 | 88.05789 | 154.54566 | 1216.2817 | 1.422397e+04 | 10.030000 | 9.77 | 155.64294 | 4.51 | 233.48458 |
| 3 | Stock total Cadmium EDTA | 1 | Bretagne | 107 | 1.052459e+01 | 17.496341 | 22.27319 | 26.98086 | 31.79325 | 36.68996 | 49.05908 | 155.2103 | 3.390700e+02 | 1.590000 | 5.61 | 56.58548 | 2.80 | 75.43579 |
| 4 | Stock total Cadmium EDTA | 1 | Centre val de Loire | 120 | 4.967508e+00 | 17.832109 | 28.20919 | 40.57518 | 48.86269 | 60.09039 | 91.32618 | 294.7506 | 1.228570e+03 | 2.910000 | 3.33 | 103.73970 | 0.83 | 123.33891 |
| 6 | Stock total Cadmium EDTA | 1 | France entiere | 1683 | 2.330283e+00 | 16.959783 | 26.39146 | 41.97102 | 57.78189 | 68.68711 | 104.80991 | 1216.2817 | 5.054430e+03 | 1.660000 | 5.76 | 129.73072 | 2.67 | 192.87670 |
| 7 | Stock total Cadmium EDTA | 1 | Grand Est | 165 | 7.676721e+00 | 26.283475 | 40.85593 | 54.38829 | 63.28832 | 76.80076 | 112.06888 | 250.5844 | 1.461820e+03 | 2.960000 | 6.06 | 124.75203 | 1.82 | 153.74899 |
| 8 | Stock total Cadmium EDTA | 1 | Haut de France | 109 | 3.654633e+01 | 58.681938 | 70.58970 | 88.07144 | 101.16829 | 110.99195 | 127.25921 | 1001.2616 | 8.715500e+03 | 7.940000 | 3.67 | 164.44249 | 1.83 | 192.87670 |
| 9 | Stock total Cadmium EDTA | 1 | Ile de France | 39 | 1.571234e+01 | 27.852365 | 36.10176 | 57.03351 | 89.68199 | 75.83248 | 103.18323 | 986.3555 | 2.475623e+04 | 24.580000 | 10.26 | 94.00138 | 5.13 | 146.36687 |
| 10 | Stock total Cadmium EDTA | 1 | Normandie | 114 | 3.935199e+00 | 22.100154 | 30.09551 | 41.79345 | 51.34070 | 62.38800 | 79.25253 | 291.9800 | 1.551470e+03 | 3.460000 | 4.39 | 103.46763 | 2.63 | 139.51148 |
| 11 | Stock total Cadmium EDTA | 1 | Nouvelle Aquitaine | 230 | 3.299297e+00 | 14.610323 | 23.86915 | 35.45740 | 54.04145 | 64.87428 | 101.54409 | 555.8309 | 3.497770e+03 | 3.740000 | 6.09 | 119.83362 | 2.61 | 168.81265 |
| 12 | Stock total Cadmium EDTA | 1 | Occitanie | 214 | 3.998951e+00 | 13.829546 | 22.83575 | 39.66342 | 65.28217 | 64.84693 | 129.05038 | 1033.9086 | 1.056294e+04 | 6.970000 | 10.28 | 127.46291 | 5.14 | 184.65885 |
| 13 | Stock total Cadmium EDTA | 1 | Pays de la Loire | 114 | 9.876922e+00 | 20.379093 | 25.36829 | 34.01998 | 39.50310 | 44.68967 | 59.29259 | 201.0140 | 5.935500e+02 | 1.920000 | 6.14 | 73.28807 | 1.75 | 98.47787 |
| 14 | Stock total Cadmium EDTA | 1 | Provence Alpes Cote D’azur | 90 | 1.031738e+01 | 22.305224 | 31.48847 | 44.49473 | 52.51896 | 63.56878 | 89.61985 | 257.9022 | 1.411630e+03 | 4.110000 | 5.56 | 99.32625 | 2.22 | 138.28833 |
| 15 | Stock total Cadmium HF | 1 | Auvergne Rhone Alpes | 222 | 1.456966e+01 | 36.448655 | 48.69528 | 75.00917 | 95.05728 | 111.13138 | 164.10016 | 491.7550 | 5.662970e+03 | 4.780000 | 7.21 | 197.84380 | 2.25 | 268.13586 |
| 16 | Stock total Cadmium HF | 1 | Bourgogne Franche Comte | 133 | 2.186465e+01 | 44.911723 | 73.69906 | 117.74304 | 187.22759 | 196.73633 | 321.81262 | 3596.3162 | 1.098298e+05 | 27.800000 | 9.02 | 343.02189 | 2.26 | 532.92192 |
| 17 | Stock total Cadmium HF | 1 | Bretagne | 107 | 1.738062e+01 | 35.302291 | 46.36642 | 54.37233 | 62.71935 | 71.08862 | 90.22443 | 288.1965 | 1.131160e+03 | 2.810000 | 5.61 | 103.20442 | 2.80 | 143.88242 |
| 18 | Stock total Cadmium HF | 1 | Centre val de Loire | 120 | 1.176579e+01 | 25.878268 | 45.91769 | 78.76878 | 99.54412 | 131.46303 | 188.24798 | 478.9697 | 5.440820e+03 | 5.550000 | 3.33 | 254.04107 | 0.83 | 316.55693 |
| 20 | Stock total Cadmium HF | 1 | France entiere | 1683 | 4.413128e+00 | 34.054113 | 52.91812 | 85.85400 | 119.14266 | 146.11333 | 221.98771 | 3596.3162 | 2.099939e+04 | 3.280000 | 5.41 | 285.35808 | 2.14 | 420.71685 |
| 21 | Stock total Cadmium HF | 1 | Grand Est | 165 | 9.222114e+00 | 56.018859 | 79.22662 | 128.79305 | 144.57980 | 182.98717 | 264.05902 | 588.1707 | 8.882740e+03 | 6.640000 | 3.64 | 319.31995 | 1.21 | 445.57027 |
| 22 | Stock total Cadmium HF | 1 | Haut de France | 109 | 6.544192e+01 | 117.067063 | 145.68592 | 167.75193 | 196.35733 | 204.43658 | 275.29481 | 1407.9960 | 1.935998e+04 | 10.510000 | 9.17 | 278.73547 | 4.59 | 369.96742 |
| 23 | Stock total Cadmium HF | 1 | Ile de France | 39 | 3.317142e+01 | 50.025185 | 77.40074 | 112.69029 | 163.54720 | 162.87837 | 257.99960 | 1327.2083 | 4.493435e+04 | 32.910000 | 10.26 | 249.66684 | 5.13 | 301.72301 |
| 24 | Stock total Cadmium HF | 1 | Normandie | 114 | 4.451820e+00 | 42.085652 | 57.06115 | 81.51849 | 94.08875 | 119.30574 | 146.06862 | 478.1885 | 4.065050e+03 | 5.770000 | 2.63 | 189.70439 | 1.75 | 269.45418 |
| 25 | Stock total Cadmium HF | 1 | Nouvelle Aquitaine | 230 | 5.770280e+00 | 28.321321 | 48.07999 | 79.02535 | 115.37820 | 144.95888 | 224.76483 | 890.9004 | 1.345744e+04 | 6.930000 | 6.96 | 287.94898 | 1.74 | 416.89143 |
| 26 | Stock total Cadmium HF | 1 | Occitanie | 214 | 4.413128e+00 | 27.851933 | 48.31998 | 86.08693 | 129.07720 | 144.22892 | 258.08699 | 1766.4990 | 3.162360e+04 | 11.510000 | 7.94 | 280.78263 | 3.27 | 404.69753 |
| 27 | Stock total Cadmium HF | 1 | Pays de la Loire | 114 | 1.944073e+01 | 33.484928 | 46.00962 | 62.63706 | 72.67208 | 84.95246 | 122.98739 | 287.5891 | 1.914690e+03 | 4.100000 | 5.26 | 137.62963 | 3.51 | 156.35807 |
| 28 | Stock total Cadmium HF | 1 | Provence Alpes Cote D’azur | 90 | 1.841466e+01 | 43.658369 | 71.77142 | 98.18161 | 113.68923 | 131.09198 | 178.84769 | 546.3982 | 5.947250e+03 | 7.960000 | 7.78 | 202.67200 | 2.22 | 285.35808 |
| 29 | Stock total Cobalt HF | 1 | Auvergne Rhone Alpes | 222 | 2.339552e+02 | 1660.649827 | 2798.04734 | 4708.48116 | 5705.63661 | 7192.99437 | 10213.59069 | 37592.4344 | 2.489524e+07 | 283.620000 | 4.95 | 13368.40200 | 2.70 | 19440.92756 |
| 30 | Stock total Cobalt HF | 1 | Bourgogne Franche Comte | 133 | 4.872628e+01 | 1780.599041 | 2978.51915 | 5894.77132 | 6552.52011 | 8524.71520 | 12078.33088 | 29755.2470 | 2.146792e+07 | 396.120000 | 2.26 | 16637.67120 | 0.75 | 21912.20220 |
| 31 | Stock total Cobalt HF | 1 | Bretagne | 107 | 3.187167e+02 | 1591.261481 | 2351.26992 | 3620.06056 | 4629.70090 | 5755.55824 | 7336.71527 | 44919.1643 | 2.169187e+07 | 474.820000 | 2.80 | 10114.77300 | 0.93 | 15211.38088 |
| 32 | Stock total Cobalt HF | 1 | Centre val de Loire | 120 | 6.877770e+01 | 1339.149905 | 2408.16278 | 4258.37088 | 4628.21177 | 6149.13450 | 7958.07927 | 20077.7049 | 9.099836e+06 | 248.910000 | 2.50 | 10975.10027 | 0.83 | 13460.41938 |
| 34 | Stock total Cobalt HF | 1 | France entiere | 1683 | 6.305694e+00 | 1413.034291 | 2549.02326 | 4559.49722 | 5245.21844 | 6936.19803 | 9456.74788 | 44919.1643 | 1.580662e+07 | 88.390000 | 3.57 | 13460.41938 | 0.89 | 20077.70494 |
| 35 | Stock total Cobalt HF | 1 | Grand Est | 165 | 2.744512e+02 | 1191.673888 | 2224.04353 | 5237.55334 | 5864.79034 | 8448.27800 | 11125.93321 | 30197.2131 | 1.942922e+07 | 296.770000 | 1.21 | 16115.82453 | 0.61 | 22521.97440 |
| 36 | Stock total Cobalt HF | 1 | Haut de France | 109 | 1.040397e+03 | 3547.522126 | 5152.07268 | 6133.47480 | 6034.11730 | 7004.04642 | 7868.44671 | 11421.9511 | 3.625833e+06 | 155.710000 | 9.17 | 9405.40088 | 0.00 | 11421.95108 |
| 37 | Stock total Cobalt HF | 1 | Ile de France | 39 | 1.396904e+03 | 2270.249911 | 3234.03759 | 3995.25627 | 4194.76288 | 5119.51359 | 6216.64315 | 9789.9348 | 2.973782e+06 | 320.080000 | 2.56 | 7748.73967 | 0.00 | 9789.93485 |
| 38 | Stock total Cobalt HF | 1 | Normandie | 114 | 6.305694e+00 | 1843.888354 | 2985.08152 | 4223.75013 | 4235.78842 | 5395.92033 | 6380.76742 | 14538.4612 | 3.999525e+06 | 185.020000 | 0.88 | 8529.37380 | 0.88 | 8529.37380 |
| 39 | Stock total Cobalt HF | 1 | Nouvelle Aquitaine | 230 | 1.093355e+02 | 826.657704 | 1621.82974 | 4374.53751 | 5281.54177 | 7420.61495 | 10434.32811 | 33869.2661 | 2.116131e+07 | 290.110000 | 3.48 | 15762.52573 | 0.43 | 22873.03290 |
| 40 | Stock total Cobalt HF | 1 | Occitanie | 214 | 3.056341e+01 | 1558.781919 | 2494.33400 | 4304.71726 | 5047.94465 | 6429.21485 | 9511.03733 | 29017.2064 | 1.305611e+07 | 231.920000 | 4.21 | 12364.36778 | 0.47 | 16303.32913 |
| 41 | Stock total Cobalt HF | 1 | Pays de la Loire | 114 | 6.472562e+02 | 1363.359325 | 2111.89875 | 3931.71375 | 4794.43564 | 6365.50071 | 9339.79089 | 18364.5629 | 1.166044e+07 | 317.780000 | 3.51 | 11045.48913 | 0.00 | 18364.56287 |
| 42 | Stock total Cobalt HF | 1 | Provence Alpes Cote D’azur | 90 | 2.974035e+02 | 1776.307722 | 2532.27775 | 3870.21041 | 4433.93169 | 6340.68301 | 8041.30306 | 10822.9177 | 6.074918e+06 | 288.810000 | 0.00 | 10822.91773 | 0.00 | 10822.91773 |
| 43 | Stock total Chrome EDTA | 1 | Auvergne Rhone Alpes | 222 | 1.243404e+00 | 14.143454 | 23.07790 | 37.44892 | 44.86518 | 59.32642 | 85.28441 | 321.5824 | 1.135240e+03 | 2.100000 | 1.80 | 114.45037 | 0.45 | 158.39212 |
| 44 | Stock total Chrome EDTA | 1 | Bourgogne Franche Comte | 133 | 2.823261e+00 | 13.238262 | 19.51248 | 30.67206 | 36.09331 | 45.29268 | 67.15768 | 121.9790 | 5.504900e+02 | 2.150000 | 6.02 | 82.56369 | 0.00 | 121.97898 |
| 45 | Stock total Chrome EDTA | 1 | Bretagne | 107 | 8.212952e+00 | 29.517638 | 43.50690 | 52.50014 | 64.68496 | 78.83839 | 103.03931 | 263.5351 | 1.512560e+03 | 3.590000 | 3.74 | 128.02333 | 1.87 | 149.88946 |
| 46 | Stock total Chrome EDTA | 1 | Centre val de Loire | 120 | 7.227164e+00 | 14.012183 | 21.19006 | 28.56923 | 37.77419 | 42.94778 | 69.59615 | 167.8057 | 7.355000e+02 | 2.330000 | 8.33 | 73.88873 | 3.33 | 102.36929 |
| 48 | Stock total Chrome EDTA | 1 | France entiere | 1683 | 6.305694e-01 | 14.016645 | 23.02388 | 36.84230 | 46.21482 | 56.97129 | 86.38377 | 514.4019 | 1.499250e+03 | 0.887603 | 5.53 | 107.57127 | 1.54 | 158.39212 |
| 49 | Stock total Chrome EDTA | 1 | Grand Est | 165 | 6.167825e+00 | 18.057745 | 25.79802 | 39.36757 | 48.32702 | 61.71725 | 87.97184 | 212.1075 | 1.127710e+03 | 2.440000 | 4.24 | 114.78325 | 1.82 | 160.78063 |
| 50 | Stock total Chrome EDTA | 1 | Haut de France | 109 | 6.751823e+00 | 18.276221 | 23.91534 | 36.67838 | 43.35034 | 50.22858 | 84.93713 | 156.2955 | 8.138000e+02 | 2.520000 | 7.34 | 89.61773 | 1.83 | 126.18133 |
| 51 | Stock total Chrome EDTA | 1 | Ile de France | 39 | 6.227750e+00 | 15.104533 | 22.46212 | 32.51243 | 53.38359 | 58.57754 | 113.59550 | 268.4357 | 2.994510e+03 | 9.090000 | 10.26 | 110.52884 | 5.13 | 133.06462 |
| 52 | Stock total Chrome EDTA | 1 | Normandie | 114 | 6.305694e-01 | 22.337249 | 35.56812 | 50.22280 | 61.07515 | 78.31622 | 106.48707 | 254.7640 | 1.825870e+03 | 3.760000 | 3.51 | 133.00236 | 2.63 | 182.14228 |
| 53 | Stock total Chrome EDTA | 1 | Nouvelle Aquitaine | 230 | 4.553129e+00 | 12.871132 | 20.44977 | 31.35922 | 39.26650 | 44.68290 | 68.98644 | 299.6770 | 1.179770e+03 | 2.290000 | 5.65 | 79.82739 | 3.04 | 110.69033 |
| 54 | Stock total Chrome EDTA | 1 | Occitanie | 214 | 7.426623e-01 | 9.464383 | 16.94018 | 30.88819 | 39.01345 | 50.27980 | 69.48587 | 471.9850 | 1.607840e+03 | 2.550000 | 2.80 | 96.33521 | 0.93 | 143.95961 |
| 55 | Stock total Chrome EDTA | 1 | Pays de la Loire | 114 | 4.903511e+00 | 15.294244 | 25.41030 | 45.47653 | 54.13686 | 69.82839 | 99.80120 | 234.5818 | 1.658170e+03 | 3.690000 | 4.39 | 133.05327 | 1.75 | 186.23278 |
| 56 | Stock total Chrome EDTA | 1 | Provence Alpes Cote D’azur | 90 | 5.424433e+00 | 15.453772 | 23.61143 | 44.38099 | 51.29236 | 70.35034 | 106.52734 | 158.3485 | 1.264030e+03 | 3.330000 | 3.33 | 124.92294 | 0.00 | 158.34855 |
| 57 | Stock total Chrome RMQS1 HF | 1 | Auvergne Rhone Alpes | 222 | 1.215306e+03 | 9503.045637 | 14081.53974 | 23067.95453 | 30109.01385 | 35422.67018 | 54406.42938 | 236463.1375 | 8.945539e+08 | 1894.720000 | 5.86 | 66501.80183 | 2.70 | 98355.86308 |
| 58 | Stock total Chrome RMQS1 HF | 1 | Bourgogne Franche Comte | 132 | 4.593395e+02 | 9907.459291 | 18183.48695 | 26952.28947 | 29784.89979 | 41050.95207 | 50324.87039 | 76243.9093 | 2.658602e+08 | 1438.770000 | 0.76 | 71176.57876 | 0.00 | 76243.90933 |
| 59 | Stock total Chrome RMQS1 HF | 1 | Bretagne | 107 | 2.446930e+03 | 11184.523786 | 16799.60286 | 25918.65867 | 35845.25805 | 31428.50014 | 38448.93988 | 1181147.4495 | 1.259250e+10 | 11048.450000 | 0.93 | 46365.60248 | 0.93 | 46365.60248 |
| 60 | Stock total Chrome RMQS1 HF | 1 | Centre val de Loire | 120 | 2.054587e+03 | 6074.120801 | 11657.63586 | 21899.92850 | 20726.99667 | 28223.82693 | 34523.61224 | 57701.6496 | 1.163326e+08 | 980.240000 | 0.83 | 46656.00785 | 0.00 | 57701.64960 |
| 62 | Stock total Chrome RMQS1 HF | 1 | France entiere | 1682 | 7.566833e+01 | 8413.474957 | 14918.48796 | 24329.35881 | 27339.46789 | 34037.91931 | 44361.48896 | 1181147.4495 | 1.181885e+09 | 812.360000 | 2.97 | 62450.02133 | 0.83 | 87586.36891 |
| 63 | Stock total Chrome RMQS1 HF | 1 | Grand Est | 165 | 1.969196e+03 | 6914.744855 | 13196.00198 | 30661.02117 | 30842.53995 | 44221.81947 | 53740.45920 | 97978.5240 | 3.753538e+08 | 1385.370000 | 1.21 | 83813.57893 | 0.00 | 97978.52400 |
| 64 | Stock total Chrome RMQS1 HF | 1 | Haut de France | 109 | 4.691949e+03 | 20805.291650 | 27194.77667 | 32360.59742 | 32411.44564 | 38130.86236 | 43528.58247 | 58959.4038 | 9.073604e+07 | 879.370000 | 4.59 | 50242.55450 | 0.00 | 58959.40376 |
| 65 | Stock total Chrome RMQS1 HF | 1 | Ile de France | 39 | 6.833129e+03 | 12553.792985 | 15889.65072 | 20563.57077 | 21550.43284 | 28011.05773 | 31978.36555 | 37117.9528 | 6.359971e+07 | 1027.700000 | 0.00 | 37117.95283 | 0.00 | 37117.95283 |
| 66 | Stock total Chrome RMQS1 HF | 1 | Normandie | 114 | 7.566833e+01 | 14065.908682 | 19461.46562 | 23378.38797 | 23586.70735 | 28653.42064 | 33109.95050 | 46880.4920 | 7.386568e+07 | 830.160000 | 6.14 | 40369.76173 | 0.00 | 46880.49205 |
| 67 | Stock total Chrome RMQS1 HF | 1 | Nouvelle Aquitaine | 230 | 1.040541e+03 | 5629.443826 | 11167.00708 | 23040.70941 | 23852.39237 | 33391.03660 | 41695.88441 | 77079.8720 | 2.245654e+08 | 943.690000 | 1.30 | 66558.12633 | 0.00 | 77079.87200 |
| 68 | Stock total Chrome RMQS1 HF | 1 | Occitanie | 214 | 1.471043e+02 | 8858.616799 | 14330.51926 | 22110.32572 | 24392.48635 | 32151.85770 | 41999.39569 | 86751.2789 | 1.926124e+08 | 878.960000 | 2.34 | 54803.07390 | 0.47 | 81270.51035 |
| 69 | Stock total Chrome RMQS1 HF | 1 | Pays de la Loire | 114 | 3.334397e+03 | 11147.586036 | 15203.04067 | 24187.02104 | 24540.23096 | 30325.08668 | 39950.94031 | 64186.0287 | 1.508323e+08 | 1156.630000 | 3.51 | 52432.91693 | 0.00 | 64186.02867 |
| 70 | Stock total Chrome RMQS1 HF | 1 | Provence Alpes Cote D’azur | 90 | 1.783524e+03 | 10037.193313 | 14654.94358 | 22735.50421 | 24091.35806 | 31138.49127 | 37364.47991 | 66524.2792 | 1.533980e+08 | 1116.930000 | 2.22 | 54766.64133 | 0.00 | 66524.27922 |
| 71 | Stock total Cuivre EDTA | 1 | Auvergne Rhone Alpes | 222 | 3.922602e+01 | 266.206284 | 406.84760 | 810.51069 | 1802.91601 | 1526.01705 | 3373.07118 | 48814.0026 | 2.390289e+07 | 323.500000 | 10.81 | 2937.67626 | 4.95 | 4774.87346 |
| 72 | Stock total Cuivre EDTA | 1 | Bourgogne Franche Comte | 133 | 3.979080e+01 | 266.088719 | 446.38355 | 945.43096 | 1391.74104 | 1537.49626 | 2397.53387 | 18688.1150 | 4.467330e+06 | 164.470000 | 6.02 | 3105.58374 | 2.26 | 3895.05584 |
| 73 | Stock total Cuivre EDTA | 1 | Bretagne | 107 | 8.278600e+01 | 310.342575 | 625.63307 | 1248.90800 | 1390.55865 | 1769.52066 | 2299.38518 | 9219.0925 | 1.452614e+06 | 100.120000 | 3.74 | 3422.03143 | 1.87 | 4306.28287 |
| 74 | Stock total Cuivre EDTA | 1 | Centre val de Loire | 120 | 7.431781e+01 | 215.679986 | 421.08693 | 709.92445 | 1082.83919 | 1176.47756 | 2132.03316 | 7793.5107 | 1.497810e+06 | 114.530000 | 8.33 | 2283.01023 | 4.17 | 3463.51277 |
| 76 | Stock total Cuivre EDTA | 1 | France entiere | 1683 | 2.749283e+01 | 285.357145 | 534.33271 | 1038.45865 | 2383.60299 | 1907.28455 | 4258.37274 | 100304.8421 | 3.281088e+07 | 132.020000 | 10.40 | 3895.05584 | 6.83 | 6026.11420 |
| 77 | Stock total Cuivre EDTA | 1 | Grand Est | 165 | 1.092941e+02 | 259.363135 | 492.65833 | 1077.28373 | 1586.47257 | 1917.86124 | 2946.44622 | 21398.8495 | 4.585299e+06 | 170.190000 | 6.06 | 3798.51656 | 2.42 | 5597.28064 |
| 78 | Stock total Cuivre EDTA | 1 | Haut de France | 109 | 1.976339e+02 | 887.254923 | 1401.07354 | 1779.03089 | 2037.14789 | 2622.99002 | 3573.56736 | 5504.0532 | 1.222115e+06 | 87.450000 | 4.59 | 4312.83804 | 0.00 | 5504.05322 |
| 79 | Stock total Cuivre EDTA | 1 | Ile de France | 39 | 3.176153e+02 | 472.144006 | 643.32698 | 1183.25394 | 3124.00483 | 2416.97708 | 6073.10325 | 36616.4747 | 3.928030e+07 | 936.350000 | 12.82 | 4892.44557 | 7.69 | 6473.24854 |
| 80 | Stock total Cuivre EDTA | 1 | Normandie | 114 | 2.749283e+01 | 406.937316 | 679.86885 | 1055.40274 | 1223.97468 | 1564.58918 | 2025.84409 | 7400.4859 | 7.937852e+05 | 54.110000 | 2.63 | 2689.74192 | 0.88 | 3607.35824 |
| 81 | Stock total Cuivre EDTA | 1 | Nouvelle Aquitaine | 230 | 7.834855e+01 | 224.128460 | 448.11041 | 895.41229 | 3978.97881 | 2047.90363 | 8784.60521 | 100304.8421 | 1.210958e+08 | 699.310000 | 14.78 | 4286.70761 | 13.04 | 6266.76532 |
| 82 | Stock total Cuivre EDTA | 1 | Occitanie | 214 | 5.127226e+01 | 299.724966 | 591.92486 | 1153.52190 | 4382.20236 | 3340.39841 | 13182.46620 | 46194.5053 | 6.118079e+07 | 420.140000 | 16.36 | 7379.25973 | 12.15 | 11444.73547 |
| 83 | Stock total Cuivre EDTA | 1 | Pays de la Loire | 114 | 2.313235e+02 | 444.341023 | 610.55987 | 934.37356 | 1772.36551 | 1632.37430 | 2604.89871 | 30953.3292 | 1.141987e+07 | 301.130000 | 7.02 | 3105.24040 | 5.26 | 3743.70289 |
| 84 | Stock total Cuivre EDTA | 1 | Provence Alpes Cote D’azur | 90 | 1.594804e+02 | 420.803598 | 678.06517 | 1541.53077 | 3497.33984 | 4601.36585 | 10318.01787 | 21927.6156 | 2.116223e+07 | 458.130000 | 10.00 | 10238.41465 | 3.33 | 14721.45316 |
| 85 | Stock total Cuivre HF | 1 | Auvergne Rhone Alpes | 222 | 3.957625e+02 | 2635.826702 | 4608.52102 | 6860.64096 | 9033.20897 | 11169.15666 | 15961.13824 | 96594.3167 | 9.037135e+07 | 606.530000 | 3.15 | 20599.56795 | 1.35 | 26590.62187 |
| 86 | Stock total Cuivre HF | 1 | Bourgogne Franche Comte | 133 | 1.262136e+02 | 2161.517279 | 3617.24173 | 6305.59920 | 7259.86412 | 9094.86400 | 12472.39075 | 39617.6110 | 3.112538e+07 | 486.080000 | 3.01 | 16532.33893 | 1.50 | 18655.93260 |
| 87 | Stock total Cuivre HF | 1 | Bretagne | 107 | 3.187167e+02 | 2929.889761 | 5069.88608 | 8155.19773 | 8750.63074 | 11460.55728 | 14752.83289 | 31088.3269 | 2.373543e+07 | 427.070000 | 0.93 | 19585.93341 | 0.93 | 19585.93341 |
| 88 | Stock total Cuivre HF | 1 | Centre val de Loire | 120 | 4.359923e+02 | 1602.878815 | 2830.34914 | 4641.65670 | 4931.59087 | 6057.43100 | 8085.50787 | 18221.9627 | 9.669414e+06 | 260.980000 | 4.17 | 10825.61581 | 0.83 | 14661.71467 |
| 90 | Stock total Cuivre HF | 1 | France entiere | 1683 | 7.175880e+01 | 2389.572942 | 4144.35490 | 6976.60886 | 9696.17037 | 11022.12266 | 17071.13634 | 216289.4696 | 1.499403e+08 | 269.850000 | 6.36 | 21272.43940 | 3.51 | 31599.65426 |
| 91 | Stock total Cuivre HF | 1 | Grand Est | 165 | 6.626844e+02 | 2053.634495 | 3666.78752 | 7659.87867 | 8303.53781 | 11027.51040 | 14522.58087 | 45493.3387 | 3.555045e+07 | 388.170000 | 1.82 | 21272.43940 | 0.61 | 29947.03774 |
| 92 | Stock total Cuivre HF | 1 | Haut de France | 109 | 1.467191e+03 | 5619.579553 | 6904.95656 | 8540.03476 | 8719.81356 | 10441.82939 | 11952.90281 | 17906.5885 | 8.303578e+06 | 253.660000 | 3.67 | 15078.34351 | 0.00 | 17906.58846 |
| 93 | Stock total Cuivre HF | 1 | Ile de France | 39 | 2.206702e+03 | 3754.666093 | 4415.58811 | 5766.94976 | 9111.73894 | 8499.71631 | 17080.12960 | 57622.4572 | 1.068759e+08 | 1724.250000 | 12.82 | 13897.61120 | 7.69 | 17721.70987 |
| 94 | Stock total Cuivre HF | 1 | Normandie | 114 | 7.175880e+01 | 2755.127508 | 4567.58679 | 5960.40689 | 6319.01567 | 7858.96251 | 9562.31297 | 19480.2983 | 9.478260e+06 | 257.060000 | 3.51 | 12555.53647 | 0.88 | 16242.72720 |
| 95 | Stock total Cuivre HF | 1 | Nouvelle Aquitaine | 230 | 3.434927e+02 | 1720.129173 | 3430.09346 | 6460.39725 | 12298.62441 | 12430.80536 | 24990.20064 | 152552.9327 | 3.734425e+08 | 1306.870000 | 10.00 | 24843.87240 | 5.65 | 37549.78673 |
| 96 | Stock total Cuivre HF | 1 | Occitanie | 214 | 1.642426e+02 | 2667.824950 | 4946.15291 | 9245.21505 | 15393.75379 | 15503.21017 | 35516.68486 | 216289.4696 | 4.336807e+08 | 1307.860000 | 11.68 | 30828.62434 | 7.48 | 46996.33707 |
| 97 | Stock total Cuivre HF | 1 | Pays de la Loire | 114 | 1.426941e+03 | 2910.506508 | 3933.21814 | 5738.23000 | 7896.53272 | 8985.88119 | 14306.00270 | 66410.7071 | 5.775787e+07 | 678.600000 | 7.02 | 15710.55178 | 4.39 | 19401.90500 |
| 98 | Stock total Cuivre HF | 1 | Provence Alpes Cote D’azur | 90 | 7.042922e+02 | 3612.471173 | 5639.89383 | 10164.60349 | 13488.90330 | 18383.85200 | 28030.14350 | 52232.0985 | 1.251934e+08 | 1276.200000 | 4.44 | 34099.65120 | 0.00 | 52232.09847 |
| 99 | Stock total Molybdene HF | 1 | Auvergne Rhone Alpes | 222 | 1.737152e+01 | 126.565345 | 204.30153 | 297.98960 | 389.99525 | 400.84606 | 654.55510 | 2869.7028 | 1.539602e+05 | 24.680000 | 9.01 | 696.11870 | 5.41 | 956.14331 |
| 100 | Stock total Molybdene HF | 1 | Bourgogne Franche Comte | 133 | 1.333341e+01 | 161.714544 | 241.33135 | 365.62715 | 468.88699 | 543.46895 | 927.84741 | 2647.3220 | 1.357816e+05 | 29.400000 | 8.27 | 992.23584 | 3.01 | 1407.25107 |
| 101 | Stock total Molybdene HF | 1 | Bretagne | 107 | 3.175649e+01 | 201.966427 | 240.75875 | 312.49020 | 348.33244 | 391.20387 | 468.99545 | 1939.4215 | 4.975078e+04 | 20.520000 | 3.74 | 611.01960 | 3.74 | 611.01960 |
| 102 | Stock total Molybdene HF | 1 | Centre val de Loire | 120 | 2.740523e+01 | 75.887799 | 133.66773 | 215.05025 | 215.07629 | 275.81596 | 343.42610 | 805.4023 | 1.371832e+04 | 9.540000 | 2.50 | 430.40289 | 0.83 | 527.15012 |
| 104 | Stock total Molybdene HF | 1 | France entiere | 1683 | 5.069778e+00 | 109.788138 | 189.99072 | 286.19639 | 355.41306 | 404.25690 | 600.09802 | 4220.3441 | 1.134374e+05 | 7.840000 | 6.65 | 725.03371 | 3.39 | 1044.72660 |
| 105 | Stock total Molybdene HF | 1 | Grand Est | 165 | 2.926660e+01 | 99.384689 | 179.75662 | 322.01282 | 450.52311 | 520.72915 | 1025.43863 | 3963.0574 | 2.390715e+05 | 36.520000 | 10.30 | 1008.35680 | 2.42 | 1493.32019 |
| 106 | Stock total Molybdene HF | 1 | Haut de France | 109 | 6.559205e+01 | 235.000287 | 293.50908 | 351.02609 | 361.77543 | 402.06110 | 489.48461 | 1221.8284 | 1.965059e+04 | 11.820000 | 10.09 | 551.10940 | 1.83 | 668.27621 |
| 107 | Stock total Molybdene HF | 1 | Ile de France | 39 | 1.246678e+02 | 157.423344 | 184.37524 | 255.01561 | 289.83267 | 379.22251 | 511.29509 | 669.0673 | 1.950199e+04 | 15.170000 | 0.00 | 669.06735 | 0.00 | 669.06735 |
| 108 | Stock total Molybdene HF | 1 | Normandie | 114 | 5.069778e+00 | 146.766858 | 195.56426 | 240.07822 | 245.53898 | 289.41197 | 350.62157 | 607.5238 | 9.025520e+03 | 7.980000 | 5.26 | 418.55757 | 0.88 | 535.74686 |
| 109 | Stock total Molybdene HF | 1 | Nouvelle Aquitaine | 230 | 2.230444e+01 | 85.274389 | 158.01858 | 282.84197 | 364.45732 | 444.00858 | 663.68366 | 4220.3441 | 1.621419e+05 | 26.210000 | 4.78 | 858.00110 | 1.74 | 1299.09289 |
| 110 | Stock total Molybdene HF | 1 | Occitanie | 214 | 5.784197e+00 | 116.529041 | 194.16791 | 284.58482 | 344.22548 | 406.94705 | 583.86743 | 1907.2464 | 7.404495e+04 | 18.010000 | 6.07 | 714.74314 | 2.80 | 1028.78906 |
| 111 | Stock total Molybdene HF | 1 | Pays de la Loire | 114 | 6.262842e+01 | 124.008578 | 196.25581 | 268.53755 | 322.21494 | 358.84550 | 512.85269 | 3110.0726 | 9.688781e+04 | 29.920000 | 6.14 | 571.57741 | 1.75 | 838.45346 |
| 112 | Stock total Molybdene HF | 1 | Provence Alpes Cote D’azur | 90 | 6.195875e+01 | 125.224958 | 180.14494 | 297.54181 | 334.36177 | 438.31465 | 609.14081 | 1192.4302 | 4.048497e+04 | 20.310000 | 1.11 | 764.33809 | 0.00 | 1192.43019 |
| 113 | Stock total Plomb EDTA | 1 | Auvergne Rhone Alpes | 222 | 8.930924e+01 | 520.049397 | 1302.14782 | 1941.91022 | 2220.99150 | 2821.35082 | 3809.97384 | 11568.7823 | 2.485017e+06 | 101.960000 | 4.05 | 5084.69568 | 1.80 | 7371.36545 |
| 114 | Stock total Plomb EDTA | 1 | Bourgogne Franche Comte | 133 | 5.100149e+02 | 1065.354555 | 1601.41397 | 2812.71028 | 3145.82673 | 3520.84267 | 4370.20742 | 35637.7620 | 1.302182e+07 | 321.270000 | 3.01 | 5361.36536 | 1.50 | 7772.29839 |
| 115 | Stock total Plomb EDTA | 1 | Bretagne | 107 | 2.917277e+02 | 922.181722 | 1107.19077 | 1401.95888 | 1676.35446 | 1802.30730 | 2306.79611 | 10844.5640 | 1.877261e+06 | 129.320000 | 3.74 | 2829.73728 | 3.74 | 2829.73728 |
| 116 | Stock total Plomb EDTA | 1 | Centre val de Loire | 120 | 5.003879e+02 | 894.792297 | 1474.39396 | 1903.18301 | 2217.17689 | 2546.34816 | 3418.31842 | 10971.0433 | 2.164810e+06 | 105.560000 | 5.83 | 3723.58350 | 2.50 | 5399.98693 |
| 118 | Stock total Plomb EDTA | 1 | France entiere | 1683 | 8.930924e+01 | 847.170416 | 1375.87830 | 2104.66250 | 2798.43848 | 3227.91939 | 4684.50695 | 57143.1758 | 1.286525e+07 | 82.590000 | 5.82 | 5985.76378 | 2.61 | 8766.02316 |
| 119 | Stock total Plomb EDTA | 1 | Grand Est | 165 | 3.029680e+02 | 731.635828 | 1410.87962 | 2575.14505 | 3240.00505 | 3369.62120 | 4605.60189 | 55301.2352 | 3.358500e+07 | 457.900000 | 4.24 | 5672.83664 | 2.42 | 8081.42580 |
| 120 | Stock total Plomb EDTA | 1 | Haut de France | 109 | 7.455417e+02 | 2013.964565 | 2437.80692 | 3361.78720 | 4386.98592 | 4534.18168 | 6879.26083 | 55637.9035 | 3.011327e+07 | 505.850000 | 8.26 | 7616.92684 | 4.59 | 10215.79970 |
| 121 | Stock total Plomb EDTA | 1 | Ile de France | 39 | 9.406402e+02 | 1163.884380 | 1721.75134 | 2610.27690 | 6860.92002 | 6603.99702 | 12867.82195 | 57143.1758 | 1.306989e+08 | 1659.370000 | 7.69 | 13295.53399 | 5.13 | 16333.78507 |
| 122 | Stock total Plomb EDTA | 1 | Normandie | 114 | 4.523123e+02 | 1165.105773 | 1650.58578 | 2000.21553 | 2524.69884 | 2617.17904 | 4095.79976 | 17510.9853 | 4.511775e+06 | 197.500000 | 10.53 | 3919.46244 | 5.26 | 4875.94351 |
| 123 | Stock total Plomb EDTA | 1 | Nouvelle Aquitaine | 230 | 1.214075e+02 | 761.620252 | 1286.53710 | 2324.68336 | 2918.81641 | 3698.48798 | 5736.44723 | 15513.5624 | 5.455555e+06 | 131.740000 | 5.65 | 7335.30689 | 1.74 | 10336.71567 |
| 124 | Stock total Plomb EDTA | 1 | Occitanie | 214 | 3.579521e+02 | 762.930855 | 1376.05865 | 2192.95121 | 2795.30225 | 3319.17577 | 5124.48525 | 20784.5291 | 6.358529e+06 | 169.790000 | 7.94 | 5756.74675 | 1.87 | 9156.96440 |
| 125 | Stock total Plomb EDTA | 1 | Pays de la Loire | 114 | 3.737923e+02 | 1050.447517 | 1438.93712 | 1880.21037 | 2276.05454 | 2541.42218 | 3633.82122 | 15499.1747 | 3.225557e+06 | 159.680000 | 7.02 | 4094.45362 | 2.63 | 5530.06612 |
| 126 | Stock total Plomb EDTA | 1 | Provence Alpes Cote D’azur | 90 | 5.109768e+02 | 828.443999 | 1138.42424 | 1877.18580 | 2513.33613 | 3193.14330 | 4697.06991 | 12959.0392 | 4.039729e+06 | 200.750000 | 4.44 | 5903.84019 | 2.22 | 7028.93420 |
| 127 | Stock total Plomb HF | 1 | Auvergne Rhone Alpes | 222 | 6.444275e+02 | 5073.285996 | 8053.96113 | 12630.69390 | 14614.02859 | 18047.87567 | 26607.88966 | 74626.3200 | 9.625469e+07 | 553.130000 | 4.05 | 32734.68213 | 0.90 | 46969.32804 |
| 128 | Stock total Plomb HF | 1 | Bourgogne Franche Comte | 133 | 8.795862e+02 | 5833.036960 | 9987.85273 | 15515.94767 | 17643.61848 | 20877.41240 | 28243.23757 | 126824.2690 | 2.368074e+08 | 1311.920000 | 2.26 | 36935.09037 | 2.26 | 36935.09037 |
| 129 | Stock total Plomb HF | 1 | Bretagne | 107 | 2.710897e+03 | 8812.438136 | 10137.15831 | 11622.88026 | 12794.75212 | 13525.80445 | 17022.81864 | 40826.0315 | 2.661832e+07 | 495.170000 | 9.35 | 18585.07285 | 3.74 | 23652.83627 |
| 130 | Stock total Plomb HF | 1 | Centre val de Loire | 120 | 1.904248e+03 | 6001.325553 | 9380.80287 | 13263.28002 | 13241.55181 | 16060.56622 | 21355.12089 | 30033.6882 | 3.148214e+07 | 510.180000 | 2.50 | 24339.77329 | 0.00 | 30033.68824 |
| 132 | Stock total Plomb HF | 1 | France entiere | 1683 | 5.637290e+02 | 5588.903738 | 8558.29977 | 12541.24717 | 14312.61213 | 17039.87555 | 23287.87521 | 214497.8450 | 1.315790e+08 | 264.660000 | 4.40 | 29593.91583 | 1.25 | 42320.26113 |
| 133 | Stock total Plomb HF | 1 | Grand Est | 165 | 2.205099e+03 | 4876.204511 | 7371.10653 | 13263.74280 | 14867.06682 | 16775.88427 | 21441.23800 | 214497.8450 | 3.992800e+08 | 1550.420000 | 1.82 | 27824.77787 | 1.21 | 38683.99871 |
| 134 | Stock total Plomb HF | 1 | Haut de France | 109 | 4.376329e+03 | 11262.083216 | 12776.97787 | 15058.99631 | 17412.64629 | 19002.07104 | 24407.35968 | 82739.4724 | 9.044411e+07 | 765.080000 | 6.42 | 26418.67312 | 2.75 | 37118.04080 |
| 135 | Stock total Plomb HF | 1 | Ile de France | 39 | 4.670427e+03 | 6037.925487 | 6645.16126 | 10086.47054 | 19788.27883 | 18079.44974 | 38299.32369 | 161809.0198 | 8.506734e+08 | 4293.080000 | 12.82 | 26546.79144 | 5.13 | 49401.66184 |
| 136 | Stock total Plomb HF | 1 | Normandie | 114 | 5.637290e+02 | 7409.265849 | 8900.08030 | 10898.63193 | 11148.93419 | 12709.62900 | 14576.63111 | 42158.3717 | 2.048613e+07 | 383.660000 | 4.39 | 17826.24507 | 1.75 | 23377.49840 |
| 137 | Stock total Plomb HF | 1 | Nouvelle Aquitaine | 230 | 1.895878e+03 | 4899.634993 | 8038.78066 | 13610.11505 | 14399.11628 | 18189.78607 | 24727.83444 | 46899.4536 | 6.929254e+07 | 489.560000 | 3.91 | 32116.39687 | 0.00 | 46899.45360 |
| 138 | Stock total Plomb HF | 1 | Occitanie | 214 | 5.755633e+02 | 5197.254927 | 8012.50401 | 12318.33635 | 13676.70664 | 16175.37103 | 22884.23590 | 66533.4129 | 7.499628e+07 | 572.040000 | 3.74 | 28332.90930 | 1.40 | 37760.95572 |
| 139 | Stock total Plomb HF | 1 | Pays de la Loire | 114 | 2.643984e+03 | 5908.595645 | 8369.97293 | 11626.52816 | 13520.11249 | 16285.08056 | 24934.43336 | 47684.1523 | 6.452162e+07 | 716.600000 | 5.26 | 28007.67763 | 1.75 | 35817.52984 |
| 140 | Stock total Plomb HF | 1 | Provence Alpes Cote D’azur | 90 | 2.708386e+03 | 5395.887867 | 7723.23395 | 10445.34996 | 12952.36122 | 16279.14923 | 22292.91810 | 41120.8561 | 6.680334e+07 | 781.260000 | 5.56 | 29085.05871 | 0.00 | 41120.85615 |
| 141 | Stock total Zinc EDTA | 1 | Auvergne Rhone Alpes | 222 | 9.897855e+01 | 311.994919 | 447.14494 | 688.88578 | 883.27983 | 971.31838 | 1785.56208 | 4738.3476 | 5.475761e+05 | 50.590000 | 10.36 | 1712.67446 | 4.95 | 2346.69372 |
| 142 | Stock total Zinc EDTA | 1 | Bourgogne Franche Comte | 133 | 1.465104e+02 | 388.972290 | 526.91808 | 798.64838 | 1405.05448 | 1102.46791 | 1405.85796 | 48632.7083 | 2.091748e+07 | 399.350000 | 6.02 | 1857.26016 | 3.01 | 2690.48893 |
| 143 | Stock total Zinc EDTA | 1 | Bretagne | 107 | 1.552522e+02 | 451.561120 | 596.50172 | 894.69420 | 1105.17920 | 1304.32034 | 2052.23449 | 4902.7364 | 5.861704e+05 | 75.190000 | 8.41 | 2184.72039 | 0.93 | 3219.72584 |
| 144 | Stock total Zinc EDTA | 1 | Centre val de Loire | 120 | 8.767335e+01 | 240.117561 | 364.17495 | 548.37268 | 632.90038 | 788.83011 | 1170.99460 | 2004.2870 | 1.546587e+05 | 36.450000 | 6.67 | 1413.54966 | 0.00 | 2004.28696 |
| 146 | Stock total Zinc EDTA | 1 | France entiere | 1683 | 8.767335e+01 | 334.361826 | 493.98854 | 756.44071 | 1198.51815 | 1219.30413 | 2083.24789 | 94353.7178 | 8.710847e+06 | 64.170000 | 8.26 | 2305.87576 | 3.74 | 3373.52440 |
| 147 | Stock total Zinc EDTA | 1 | Grand Est | 165 | 2.745379e+02 | 537.091276 | 679.26605 | 904.48867 | 1214.36477 | 1346.24851 | 2036.05924 | 12996.3580 | 1.470733e+06 | 97.150000 | 6.67 | 2235.09972 | 3.03 | 3287.77742 |
| 148 | Stock total Zinc EDTA | 1 | Haut de France | 109 | 4.901646e+02 | 855.500734 | 1171.20222 | 1650.81224 | 2188.64913 | 2454.94455 | 3380.17150 | 23725.5039 | 6.109272e+06 | 210.550000 | 4.59 | 4126.04845 | 2.75 | 5617.76560 |
| 149 | Stock total Zinc EDTA | 1 | Ile de France | 39 | 2.937915e+02 | 399.654614 | 517.87667 | 959.51941 | 4293.02629 | 1885.19648 | 4622.40145 | 94353.7178 | 2.274842e+08 | 2391.130000 | 12.82 | 3135.85994 | 7.69 | 5629.60408 |
| 150 | Stock total Zinc EDTA | 1 | Normandie | 114 | 2.280398e+02 | 460.291998 | 653.57847 | 872.49702 | 1205.94529 | 1276.78011 | 1791.06282 | 12187.1898 | 1.864076e+06 | 126.230000 | 7.89 | 1883.84438 | 4.39 | 2860.19648 |
| 151 | Stock total Zinc EDTA | 1 | Nouvelle Aquitaine | 230 | 1.140785e+02 | 306.213210 | 396.36799 | 583.72941 | 925.88360 | 1013.80241 | 1715.95169 | 11596.2033 | 1.389191e+06 | 76.020000 | 8.70 | 1865.15843 | 3.48 | 2778.49579 |
| 152 | Stock total Zinc EDTA | 1 | Occitanie | 214 | 1.116391e+02 | 291.747314 | 411.68787 | 595.47451 | 1112.42184 | 966.20162 | 1900.84598 | 20047.3156 | 3.718529e+06 | 134.440000 | 10.75 | 1740.52684 | 7.94 | 2536.51846 |
| 153 | Stock total Zinc EDTA | 1 | Pays de la Loire | 114 | 2.234614e+02 | 370.434771 | 523.72156 | 763.10336 | 1031.48646 | 1167.03926 | 1882.85522 | 8054.5490 | 9.613739e+05 | 94.470000 | 7.02 | 2112.19384 | 2.63 | 2865.86955 |
| 154 | Stock total Zinc EDTA | 1 | Provence Alpes Cote D’azur | 90 | 1.787957e+02 | 418.793388 | 584.28452 | 865.24160 | 1255.40434 | 1368.74976 | 2406.04709 | 7494.0771 | 1.456951e+06 | 136.410000 | 8.89 | 2542.45404 | 5.56 | 3602.19912 |
| 155 | Stock total Zinc HF | 1 | Auvergne Rhone Alpes | 222 | 2.227478e+03 | 16188.869828 | 23120.22994 | 32074.21336 | 35627.70564 | 43251.26189 | 56453.46380 | 279695.2383 | 5.273654e+08 | 1579.250000 | 3.60 | 72303.36933 | 0.45 | 89317.16403 |
| 156 | Stock total Zinc HF | 1 | Bourgogne Franche Comte | 133 | 1.940674e+03 | 15438.123365 | 22764.04795 | 36755.58356 | 44840.00913 | 50642.30400 | 70186.98856 | 646923.1888 | 3.403141e+09 | 5059.450000 | 4.51 | 90032.00800 | 1.50 | 121178.72000 |
| 157 | Stock total Zinc HF | 1 | Bretagne | 107 | 5.146211e+03 | 19752.900773 | 24375.09053 | 31138.42717 | 32598.12520 | 40055.14210 | 47170.75255 | 74733.9840 | 1.335309e+08 | 1027.740000 | 0.93 | 60813.60400 | 0.00 | 74733.98400 |
| 158 | Stock total Zinc HF | 1 | Centre val de Loire | 120 | 2.440726e+03 | 9067.588701 | 13341.67128 | 20227.22414 | 22524.43628 | 29667.82495 | 37090.63748 | 64046.8960 | 1.500560e+08 | 1034.920000 | 2.50 | 49327.00600 | 0.00 | 64046.89602 |
| 160 | Stock total Zinc HF | 1 | France entiere | 1683 | 6.923652e+02 | 13256.711183 | 20433.74352 | 30505.31020 | 34502.71504 | 42029.99226 | 56533.63622 | 646923.1888 | 8.762037e+08 | 706.130000 | 3.86 | 74212.28733 | 0.83 | 104208.19813 |
| 161 | Stock total Zinc HF | 1 | Grand Est | 165 | 3.505532e+03 | 13332.083880 | 21741.36420 | 35344.35480 | 37971.28926 | 49731.24533 | 65841.47164 | 170243.2733 | 4.961922e+08 | 1485.220000 | 1.21 | 87866.81000 | 0.61 | 93355.33582 |
| 162 | Stock total Zinc HF | 1 | Haut de France | 109 | 1.323788e+04 | 28447.833073 | 32444.52000 | 38634.23240 | 39091.84398 | 44252.11031 | 50157.28021 | 119021.6928 | 1.606480e+08 | 970.960000 | 5.50 | 57784.81973 | 0.92 | 72067.06970 |
| 163 | Stock total Zinc HF | 1 | Ile de France | 39 | 7.087372e+03 | 16591.032400 | 18609.53908 | 25679.90842 | 34329.05454 | 35680.71589 | 50160.93122 | 259336.0927 | 1.581099e+09 | 5856.820000 | 7.69 | 56508.90000 | 2.56 | 73760.29645 |
| 164 | Stock total Zinc HF | 1 | Normandie | 114 | 6.923652e+02 | 15140.483557 | 21750.57197 | 26801.20193 | 26871.77109 | 32120.28594 | 38572.50556 | 55745.2992 | 1.059979e+08 | 703.750000 | 8.77 | 45338.93371 | 0.00 | 55745.29917 |
| 165 | Stock total Zinc HF | 1 | Nouvelle Aquitaine | 230 | 1.649648e+03 | 9822.034108 | 17097.79067 | 28711.02457 | 32417.16305 | 43697.62333 | 60260.30462 | 104208.1981 | 4.340071e+08 | 1295.430000 | 3.04 | 82094.92900 | 0.00 | 104208.19813 |
| 166 | Stock total Zinc HF | 1 | Occitanie | 214 | 8.254977e+02 | 14073.048129 | 19615.02570 | 32613.47206 | 38636.82764 | 44929.05807 | 65019.03812 | 623538.2992 | 2.199202e+09 | 3133.100000 | 3.74 | 79737.90242 | 1.40 | 118384.89067 |
| 167 | Stock total Zinc HF | 1 | Pays de la Loire | 114 | 6.559999e+03 | 15398.384026 | 19894.47090 | 26204.81841 | 31288.20496 | 40663.33628 | 52156.71846 | 96111.4240 | 2.994844e+08 | 1543.350000 | 3.51 | 71036.97675 | 0.00 | 96111.42401 |
| 168 | Stock total Zinc HF | 1 | Provence Alpes Cote D’azur | 90 | 4.528028e+03 | 15495.899609 | 20953.32628 | 30342.48216 | 35310.75086 | 46492.11852 | 56547.15169 | 179168.0995 | 5.080463e+08 | 2324.520000 | 2.22 | 66642.37747 | 1.11 | 88874.42400 |